{
  "cells": [
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "view-in-github",
        "colab_type": "text"
      },
      "source": [
        "<a href=\"https://colab.research.google.com/github/vitaldb/examples/blob/master/ci_aki.ipynb\" target=\"_parent\"><img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/></a>"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "bpQy0pzwVUBL"
      },
      "source": [
        "# Calculation of AKI risk depending on BP during surgery\n",
        "In this example, we will calculate the risk of acute kidney injury (AKI) depending on blood pressure during surgery, using the VitalDB open dataset."
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "SqQZp4YdV2tg"
      },
      "source": [
        "## Required libraries and datasets"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 1,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "9wmjfMpGN6_z",
        "outputId": "02efb940-90e9-4c74-cf42-a7604438ee88"
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "Looking in indexes: https://pypi.org/simple, https://us-python.pkg.dev/colab-wheels/public/simple/\n",
            "Collecting vitaldb\n",
            "  Downloading vitaldb-1.4.2-py3-none-any.whl (54 kB)\n",
            "\u001b[2K     \u001b[90m━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━\u001b[0m \u001b[32m55.0/55.0 kB\u001b[0m \u001b[31m2.7 MB/s\u001b[0m eta \u001b[36m0:00:00\u001b[0m\n",
            "\u001b[?25hRequirement already satisfied: numpy in /usr/local/lib/python3.10/dist-packages (from vitaldb) (1.22.4)\n",
            "Requirement already satisfied: pandas in /usr/local/lib/python3.10/dist-packages (from vitaldb) (1.5.3)\n",
            "Requirement already satisfied: requests in /usr/local/lib/python3.10/dist-packages (from vitaldb) (2.27.1)\n",
            "Collecting wfdb (from vitaldb)\n",
            "  Downloading wfdb-4.1.2-py3-none-any.whl (159 kB)\n",
            "\u001b[2K     \u001b[90m━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━\u001b[0m \u001b[32m160.0/160.0 kB\u001b[0m \u001b[31m6.6 MB/s\u001b[0m eta \u001b[36m0:00:00\u001b[0m\n",
            "\u001b[?25hRequirement already satisfied: python-dateutil>=2.8.1 in /usr/local/lib/python3.10/dist-packages (from pandas->vitaldb) (2.8.2)\n",
            "Requirement already satisfied: pytz>=2020.1 in /usr/local/lib/python3.10/dist-packages (from pandas->vitaldb) (2022.7.1)\n",
            "Requirement already satisfied: urllib3<1.27,>=1.21.1 in /usr/local/lib/python3.10/dist-packages (from requests->vitaldb) (1.26.15)\n",
            "Requirement already satisfied: certifi>=2017.4.17 in /usr/local/lib/python3.10/dist-packages (from requests->vitaldb) (2022.12.7)\n",
            "Requirement already satisfied: charset-normalizer~=2.0.0 in /usr/local/lib/python3.10/dist-packages (from requests->vitaldb) (2.0.12)\n",
            "Requirement already satisfied: idna<4,>=2.5 in /usr/local/lib/python3.10/dist-packages (from requests->vitaldb) (3.4)\n",
            "Requirement already satisfied: SoundFile>=0.10.0 in /usr/local/lib/python3.10/dist-packages (from wfdb->vitaldb) (0.12.1)\n",
            "Requirement already satisfied: matplotlib>=3.2.2 in /usr/local/lib/python3.10/dist-packages (from wfdb->vitaldb) (3.7.1)\n",
            "Requirement already satisfied: scipy>=1.0.0 in /usr/local/lib/python3.10/dist-packages (from wfdb->vitaldb) (1.10.1)\n",
            "Requirement already satisfied: contourpy>=1.0.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=3.2.2->wfdb->vitaldb) (1.0.7)\n",
            "Requirement already satisfied: cycler>=0.10 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=3.2.2->wfdb->vitaldb) (0.11.0)\n",
            "Requirement already satisfied: fonttools>=4.22.0 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=3.2.2->wfdb->vitaldb) (4.39.3)\n",
            "Requirement already satisfied: kiwisolver>=1.0.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=3.2.2->wfdb->vitaldb) (1.4.4)\n",
            "Requirement already satisfied: packaging>=20.0 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=3.2.2->wfdb->vitaldb) (23.1)\n",
            "Requirement already satisfied: pillow>=6.2.0 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=3.2.2->wfdb->vitaldb) (8.4.0)\n",
            "Requirement already satisfied: pyparsing>=2.3.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=3.2.2->wfdb->vitaldb) (3.0.9)\n",
            "Requirement already satisfied: six>=1.5 in /usr/local/lib/python3.10/dist-packages (from python-dateutil>=2.8.1->pandas->vitaldb) (1.16.0)\n",
            "Requirement already satisfied: cffi>=1.0 in /usr/local/lib/python3.10/dist-packages (from SoundFile>=0.10.0->wfdb->vitaldb) (1.15.1)\n",
            "Requirement already satisfied: pycparser in /usr/local/lib/python3.10/dist-packages (from cffi>=1.0->SoundFile>=0.10.0->wfdb->vitaldb) (2.21)\n",
            "Installing collected packages: wfdb, vitaldb\n",
            "Successfully installed vitaldb-1.4.2 wfdb-4.1.2\n"
          ]
        }
      ],
      "source": [
        "!pip install vitaldb\n",
        "import vitaldb\n",
        "import pandas as pd\n",
        "import numpy as np\n",
        "\n",
        "df_cases = pd.read_csv(\"https://api.vitaldb.net/cases\")  # Load clinical data\n",
        "df_trks = pd.read_csv('https://api.vitaldb.net/trks')  # Load track list\n",
        "df_labs = pd.read_csv('https://api.vitaldb.net/labs')  # Load lab result"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "lbUrsNBzV8ok"
      },
      "source": [
        "## Case Selection\n",
        "- Let's use only 100 emergency surgical patients."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 16,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "WXeQT1UPWER5",
        "outputId": "c5fef9da-1e39-4a06-d4ba-1eabfd7e3dbf"
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "Total 981 cases found\n"
          ]
        }
      ],
      "source": [
        "# inclusion / exclusion criteria\n",
        "caseids = list(\n",
        "    set(df_trks.loc[df_trks['tname'] == 'Vigileo/CI', 'caseid']) |\n",
        "    set(df_trks.loc[df_trks['tname'] == 'EV1000/CI', 'caseid']) |\n",
        "    set(df_trks.loc[df_trks['tname'] == 'Vigilance/CI', 'caseid'])\n",
        ")\n",
        "print('Total {} cases found'.format(len(caseids)))"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "OkpKkCKbWL9B"
      },
      "source": [
        "## Load and calculate data for each case\n",
        "- First of all, calculate the preoperative and postoperative creatinine concentration for each case, and then calculate the possibility that AKI occurs.\n",
        "\n",
        "- Secondly, obtain the ART_MBP data from vitaldb datasets. Then, calculate the proportion of surgeries that have a measured value lower than the threshold between 40-80 mmHg in all surgeries."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 20,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "yckBVdfxWQEz",
        "outputId": "bda4e4ad-0ac0-4775-8541-7d471dac8496",
        "scrolled": true
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "loading 4098...0.84 -> 1.14, no AKI, min CI=0.10000000149011612, max CI=3.9000000953674316\n",
            "loading 6147...no preop cr\n",
            "loading 2057...1.23 -> 1.35, no AKI, min CI=1.2999999523162842, max CI=3.200000047683716\n",
            "loading 12...0.58 -> 0.65, no AKI, min CI=3.299999952316284, max CI=3.9000000953674316\n",
            "loading 6160...no postop cr\n",
            "loading 17...1.05 -> 1.09, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 4112...no postop cr\n",
            "loading 19...1.19 -> 0.78, no AKI, min CI=0.10000000149011612, max CI=3.9000000953674316\n",
            "loading 6163...0.77 -> 0.76, no AKI, min CI=0.4000000059604645, max CI=3.9000000953674316\n",
            "loading 22...no preop cr\n",
            "loading 4120...1.03 -> 0.96, no AKI, min CI=2.9000000953674316, max CI=3.9000000953674316\n",
            "loading 2072...no postop cr\n",
            "loading 26...0.73 -> 0.78, no AKI, min CI=1.5, max CI=3.5999999046325684\n",
            "loading 25...no preop cr\n",
            "loading 29...0.78 -> 0.74, no AKI, min CI=1.0, max CI=3.0999999046325684\n",
            "loading 2082...1.46 -> 1.29, no AKI, min CI=2.0999999046325684, max CI=3.700000047683716\n",
            "loading 34...no postop cr\n",
            "loading 6180...1.31 -> 1.14, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 4133...no preop cr\n",
            "loading 38...0.95 -> 0.99, no AKI, min CI=2.4000000953674316, max CI=3.9000000953674316\n",
            "loading 6182...no preop cr\n",
            "loading 6185...0.55 -> 0.54, no AKI, min CI=2.0, max CI=3.5\n",
            "loading 6186...0.74 -> 0.64, no AKI, min CI=0.800000011920929, max CI=3.0\n",
            "loading 4140...0.89 -> 0.94, no AKI, min CI=1.0, max CI=3.700000047683716\n",
            "loading 4143...0.85 -> 1.18, no AKI, min CI=0.699999988079071, max CI=3.9000000953674316\n",
            "loading 6194...1.1 -> 1.13, no AKI, min CI=1.0, max CI=3.4000000953674316\n",
            "loading 52...6.85 -> 5.95, no AKI, min CI=2.5999999046325684, max CI=3.9000000953674316\n",
            "loading 4149...1.07 -> 1.05, no AKI, min CI=1.2999999523162842, max CI=3.5\n",
            "loading 6196...0.82 -> 1.0, no AKI, min CI=0.699999988079071, max CI=2.799999952316284\n",
            "loading 6199...0.66 -> 0.67, no AKI, min CI=2.0, max CI=3.4000000953674316\n",
            "loading 6200...1.0 -> 0.85, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 4150...no postop cr\n",
            "loading 55...0.43 -> 0.63, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 2106...0.94 -> 0.88, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 60...9.36 -> 7.95, no AKI, min CI=2.5, max CI=3.9000000953674316\n",
            "loading 2114...0.7 -> 0.61, no AKI, min CI=1.7999999523162842, max CI=2.700000047683716\n",
            "loading 6210...0.77 -> 0.73, no AKI, min CI=2.4000000953674316, max CI=3.9000000953674316\n",
            "loading 68...1.01 -> 1.11, no AKI, min CI=1.7999999523162842, max CI=3.0999999046325684\n",
            "loading 70...1.12 -> 1.05, no AKI, min CI=1.600000023841858, max CI=3.5\n",
            "loading 2121...0.7 -> 0.72, no AKI, min CI=3.0, max CI=3.9000000953674316\n",
            "loading 76...no preop cr\n",
            "loading 4172...no postop cr\n",
            "loading 6220...0.68 -> 2.33, AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 4175...0.76 -> 0.73, no AKI, min CI=2.700000047683716, max CI=3.9000000953674316\n",
            "loading 2128...0.75 -> 0.68, no AKI, min CI=1.100000023841858, max CI=2.700000047683716\n",
            "loading 83...8.13 -> 7.36, no AKI, min CI=0.699999988079071, max CI=3.200000047683716\n",
            "loading 4179...1.0 -> 1.09, no AKI, min CI=2.5999999046325684, max CI=3.9000000953674316\n",
            "loading 6227...0.76 -> 0.9, no AKI, min CI=0.5, max CI=3.9000000953674316\n",
            "loading 6228...0.64 -> 0.75, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 4183...0.72 -> 0.76, no AKI, min CI=2.9000000953674316, max CI=3.799999952316284\n",
            "loading 4187...0.93 -> 1.04, no AKI, min CI=1.2999999523162842, max CI=3.9000000953674316\n",
            "loading 6235...0.8 -> 0.74, no AKI, min CI=1.2000000476837158, max CI=3.4000000953674316\n",
            "loading 2139...no postop cr\n",
            "loading 6238...0.95 -> 0.96, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 94...0.87 -> 1.08, no AKI, min CI=1.2000000476837158, max CI=3.200000047683716\n",
            "loading 97...8.63 -> 7.0, no AKI, min CI=2.5, max CI=3.9000000953674316\n",
            "loading 6246...0.39 -> 0.42, no AKI, min CI=3.200000047683716, max CI=3.9000000953674316\n",
            "loading 104...0.92 -> 0.82, no AKI, min CI=2.4000000953674316, max CI=3.9000000953674316\n",
            "loading 4201...0.82 -> 1.03, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 4202...0.61 -> 0.72, no AKI, min CI=3.299999952316284, max CI=3.9000000953674316\n",
            "loading 2157...0.78 -> 0.74, no AKI, min CI=2.5999999046325684, max CI=3.9000000953674316\n",
            "loading 111...13.17 -> 11.77, no AKI, min CI=0.8999999761581421, max CI=3.9000000953674316\n",
            "loading 4207...0.65 -> 0.62, no AKI, min CI=1.7999999523162842, max CI=2.9000000953674316\n",
            "loading 2161...0.69 -> 0.67, no AKI, min CI=1.2999999523162842, max CI=3.9000000953674316\n",
            "loading 6257...0.52 -> 0.79, AKI, min CI=0.699999988079071, max CI=2.5\n",
            "loading 118...0.95 -> 1.02, no AKI, min CI=0.800000011920929, max CI=3.799999952316284\n",
            "loading 4214...0.77 -> 0.69, no AKI, min CI=1.100000023841858, max CI=2.4000000953674316\n",
            "loading 2168...5.43 -> 3.5, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 6262...no preop cr\n",
            "loading 6264...1.21 -> 1.47, no AKI, min CI=1.600000023841858, max CI=3.799999952316284\n",
            "loading 125...1.18 -> 1.31, no AKI, min CI=0.10000000149011612, max CI=3.9000000953674316\n",
            "loading 2174...5.42 -> 4.63, no AKI, min CI=1.5, max CI=3.5\n",
            "loading 130...0.92 -> 1.08, no AKI, min CI=1.7000000476837158, max CI=3.5999999046325684\n",
            "loading 6275...1.63 -> 1.6, no AKI, min CI=1.399999976158142, max CI=3.0\n",
            "loading 6277...0.57 -> 0.54, no AKI, min CI=2.799999952316284, max CI=3.9000000953674316\n",
            "loading 136...0.9 -> 0.82, no AKI, min CI=1.2999999523162842, max CI=3.0\n",
            "loading 6281...0.72 -> 0.77, no AKI, min CI=0.8999999761581421, max CI=3.9000000953674316\n",
            "loading 6280...0.75 -> 0.77, no AKI, min CI=0.10000000149011612, max CI=3.9000000953674316\n",
            "loading 137...no preop cr\n",
            "loading 138...0.72 -> 0.65, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 4233...13.67 -> 11.49, no AKI, min CI=1.2999999523162842, max CI=3.9000000953674316\n",
            "loading 4236...no preop cr\n",
            "loading 143...no preop cr\n",
            "loading 2185...1.03 -> 0.99, no AKI, min CI=3.700000047683716, max CI=3.799999952316284\n",
            "loading 145...2.36 -> 2.36, no AKI, min CI=3.5, max CI=3.9000000953674316\n",
            "loading 4242...no preop cr\n",
            "loading 4245...0.71 -> 0.71, no AKI, min CI=1.2999999523162842, max CI=2.9000000953674316\n",
            "loading 6293...0.94 -> 1.09, no AKI, min CI=0.30000001192092896, max CI=3.799999952316284\n",
            "loading 6295...0.63 -> 0.55, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 6296...1.05 -> 1.23, no AKI, min CI=1.100000023841858, max CI=3.9000000953674316\n",
            "loading 153...0.89 -> 0.69, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 2202...7.12 -> 3.79, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 6297...0.79 -> 0.63, no AKI, min CI=2.799999952316284, max CI=3.9000000953674316\n",
            "loading 156...0.56 -> 0.48, no AKI, min CI=1.5, max CI=3.0999999046325684\n",
            "loading 157...0.59 -> 0.51, no AKI, min CI=1.7999999523162842, max CI=3.0999999046325684\n",
            "loading 158...0.61 -> 0.68, no AKI, min CI=0.800000011920929, max CI=2.5999999046325684\n",
            "loading 4252...0.74 -> 0.73, no AKI, min CI=1.5, max CI=3.200000047683716\n",
            "loading 4255...0.45 -> 0.87, AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 6305...0.95 -> 0.94, no AKI, min CI=1.2000000476837158, max CI=3.0\n",
            "loading 4258...0.57 -> 0.56, no AKI, min CI=2.4000000953674316, max CI=3.9000000953674316\n",
            "loading 4259...0.6 -> 0.5, no AKI, min CI=2.200000047683716, max CI=3.9000000953674316\n",
            "loading 164...no ci\n",
            "loading 6309...0.78 -> 1.2, AKI, min CI=1.2000000476837158, max CI=3.200000047683716\n",
            "loading 4265...no preop cr\n",
            "loading 2221...0.98 -> 1.01, no AKI, min CI=3.0999999046325684, max CI=3.9000000953674316\n",
            "loading 2222...0.55 -> 0.51, no AKI, min CI=3.299999952316284, max CI=3.9000000953674316\n",
            "loading 175...no preop cr\n",
            "loading 2224...no postop cr\n",
            "loading 2225...no preop cr\n",
            "loading 180...no postop cr\n",
            "loading 2229...0.79 -> 0.78, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 2231...1.17 -> 0.87, no AKI, min CI=1.600000023841858, max CI=3.0999999046325684\n",
            "loading 4279...0.73 -> 0.85, no AKI, min CI=2.200000047683716, max CI=3.0999999046325684\n",
            "loading 4281...0.83 -> 0.98, no AKI, min CI=2.4000000953674316, max CI=3.9000000953674316\n",
            "loading 186...0.68 -> 0.7, no AKI, min CI=2.0999999046325684, max CI=3.200000047683716\n",
            "loading 6330...1.0 -> 0.93, no AKI, min CI=0.8999999761581421, max CI=3.9000000953674316\n",
            "loading 4283...0.75 -> 1.03, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 189...0.57 -> 0.54, no AKI, min CI=0.699999988079071, max CI=3.299999952316284\n",
            "loading 4287...0.68 -> 0.81, no AKI, min CI=1.899999976158142, max CI=3.5999999046325684\n",
            "loading 4289...0.8 -> 0.65, no AKI, min CI=1.399999976158142, max CI=3.9000000953674316\n",
            "loading 195...4.07 -> 3.3, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 4294...0.95 -> 1.03, no AKI, min CI=1.0, max CI=3.299999952316284\n",
            "loading 2248...0.65 -> 0.52, no AKI, min CI=1.2999999523162842, max CI=2.700000047683716\n",
            "loading 2249...1.28 -> 1.17, no AKI, min CI=1.399999976158142, max CI=2.9000000953674316\n",
            "loading 6345...no preop cr\n",
            "loading 202...0.88 -> 0.85, no AKI, min CI=2.9000000953674316, max CI=3.9000000953674316\n",
            "loading 4299...5.21 -> 4.46, no AKI, min CI=1.399999976158142, max CI=3.799999952316284\n",
            "loading 6346...9.8 -> 7.58, no AKI, min CI=2.5999999046325684, max CI=3.700000047683716\n",
            "loading 4302...0.66 -> 0.74, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 6351...0.75 -> 0.81, no AKI, min CI=2.200000047683716, max CI=3.9000000953674316\n",
            "loading 2252...1.04 -> 1.07, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 4305...4.94 -> 1.72, no AKI, min CI=2.5999999046325684, max CI=3.9000000953674316\n",
            "loading 210...no preop cr\n",
            "loading 6355...1.12 -> 0.93, no AKI, min CI=1.399999976158142, max CI=3.9000000953674316\n",
            "loading 4310...0.71 -> 0.84, no AKI, min CI=0.699999988079071, max CI=3.9000000953674316\n",
            "loading 217...no preop cr\n",
            "loading 218...0.87 -> 0.93, no AKI, min CI=2.5999999046325684, max CI=3.9000000953674316\n",
            "loading 2265...0.84 -> 0.66, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 2267...0.61 -> 0.79, no AKI, min CI=0.5, max CI=3.9000000953674316\n",
            "loading 2272...1.2 -> 1.79, no AKI, min CI=0.4000000059604645, max CI=3.9000000953674316\n",
            "loading 2273...0.95 -> 0.94, no AKI, min CI=1.600000023841858, max CI=3.9000000953674316\n",
            "loading 234...0.89 -> 0.75, no AKI, min CI=1.7999999523162842, max CI=3.5\n",
            "loading 2282...0.75 -> 0.76, no AKI, min CI=2.700000047683716, max CI=3.9000000953674316\n",
            "loading 236...7.05 -> 6.54, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 237...0.87 -> 0.72, no AKI, min CI=0.10000000149011612, max CI=3.9000000953674316\n",
            "loading 4333...0.71 -> 0.6, no AKI, min CI=1.0, max CI=2.799999952316284\n",
            "loading 6378...0.61 -> 0.69, no AKI, min CI=0.8999999761581421, max CI=2.9000000953674316\n",
            "loading 6383...6.7 -> 6.14, no AKI, min CI=3.0999999046325684, max CI=3.9000000953674316\n",
            "loading 241...0.67 -> 0.85, no AKI, min CI=0.30000001192092896, max CI=3.700000047683716\n",
            "loading 4341...no preop cr\n",
            "loading 251...1.01 -> 1.16, no AKI, min CI=0.5, max CI=3.9000000953674316\n",
            "loading 4350...0.97 -> 1.03, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 4352...0.62 -> 0.59, no AKI, min CI=2.0999999046325684, max CI=3.299999952316284\n",
            "loading 2304...0.75 -> 0.76, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 2307...1.01 -> 0.85, no AKI, min CI=1.2999999523162842, max CI=3.4000000953674316\n",
            "loading 4356...0.65 -> 0.57, no AKI, min CI=2.0999999046325684, max CI=3.799999952316284\n",
            "loading 2309...0.58 -> 0.58, no AKI, min CI=1.600000023841858, max CI=3.0\n",
            "loading 2311...0.73 -> 0.73, no AKI, min CI=1.399999976158142, max CI=2.700000047683716\n",
            "loading 263...0.72 -> 0.84, no AKI, min CI=1.2999999523162842, max CI=3.0999999046325684\n",
            "loading 2317...0.63 -> 0.78, no AKI, min CI=1.2000000476837158, max CI=3.200000047683716\n",
            "loading 269...no preop cr\n",
            "loading 4371...0.94 -> 1.13, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 2325...6.63 -> 4.92, no AKI, min CI=1.5, max CI=3.5\n",
            "loading 2326...0.7 -> 0.7, no AKI, min CI=2.9000000953674316, max CI=3.9000000953674316\n",
            "loading 2327...0.61 -> 1.0, AKI, min CI=2.200000047683716, max CI=3.9000000953674316\n",
            "loading 282...0.88 -> 0.86, no AKI, min CI=3.0, max CI=3.9000000953674316\n",
            "loading 284...6.82 -> 5.98, no AKI, min CI=1.5, max CI=3.4000000953674316\n",
            "loading 2332...0.57 -> 0.61, no AKI, min CI=1.399999976158142, max CI=3.9000000953674316\n",
            "loading 286...0.83 -> 1.0, no AKI, min CI=0.8999999761581421, max CI=3.9000000953674316\n",
            "loading 2335...0.58 -> 0.59, no AKI, min CI=0.699999988079071, max CI=3.700000047683716\n",
            "loading 2336...no preop cr\n",
            "loading 4383...0.67 -> 0.63, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 4386...no ci\n",
            "loading 4388...1.26 -> 1.25, no AKI, min CI=0.5, max CI=2.700000047683716\n",
            "loading 293...0.88 -> 0.83, no AKI, min CI=1.0, max CI=3.9000000953674316\n",
            "loading 4389...0.71 -> 0.67, no AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 296...0.82 -> 0.73, no AKI, min CI=2.4000000953674316, max CI=3.799999952316284\n",
            "loading 2344...0.8 -> 0.85, no AKI, min CI=0.800000011920929, max CI=3.9000000953674316\n",
            "loading 4398...no postop cr\n",
            "loading 4400...0.98 -> 0.91, no AKI, min CI=1.600000023841858, max CI=3.9000000953674316\n",
            "loading 2352...no preop cr\n",
            "loading 306...0.84 -> 0.85, no AKI, min CI=1.2999999523162842, max CI=3.9000000953674316\n",
            "loading 2353...1.19 -> 1.74, no AKI, min CI=1.600000023841858, max CI=3.200000047683716\n",
            "loading 308...0.78 -> 0.82, no AKI, min CI=0.800000011920929, max CI=2.799999952316284\n",
            "loading 2357...0.67 -> 0.5, no AKI, min CI=1.5, max CI=2.9000000953674316\n",
            "loading 4405...0.95 -> 1.0, no AKI, min CI=2.200000047683716, max CI=3.9000000953674316\n",
            "loading 2359...0.75 -> 0.51, no AKI, min CI=2.5, max CI=3.9000000953674316\n",
            "loading 2356...0.77 -> 0.88, no AKI, min CI=1.399999976158142, max CI=3.5999999046325684\n",
            "loading 4409...1.01 -> 1.03, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 318...0.87 -> 0.82, no AKI, min CI=2.9000000953674316, max CI=3.9000000953674316\n",
            "loading 321...1.28 -> 1.16, no AKI, min CI=2.5, max CI=3.9000000953674316\n",
            "loading 2372...0.95 -> 1.03, no AKI, min CI=2.4000000953674316, max CI=3.9000000953674316\n",
            "loading 2375...0.65 -> 0.58, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 4425...0.87 -> 1.07, no AKI, min CI=2.700000047683716, max CI=3.9000000953674316\n",
            "loading 2380...no postop cr\n",
            "loading 342...0.85 -> 0.88, no AKI, min CI=1.0, max CI=3.9000000953674316\n",
            "loading 343...0.54 -> 0.62, no AKI, min CI=1.2999999523162842, max CI=3.5\n",
            "loading 2392...0.74 -> 0.62, no AKI, min CI=2.799999952316284, max CI=3.9000000953674316\n",
            "loading 345...0.74 -> 0.69, no AKI, min CI=1.5, max CI=3.4000000953674316\n",
            "loading 349...1.07 -> 1.61, AKI, min CI=1.2000000476837158, max CI=3.9000000953674316\n",
            "loading 4449...0.21 -> 0.3, no AKI, min CI=2.5999999046325684, max CI=3.9000000953674316\n",
            "loading 358...0.66 -> 0.72, no AKI, min CI=1.600000023841858, max CI=2.5999999046325684\n",
            "loading 363...7.38 -> 5.97, no AKI, min CI=1.7000000476837158, max CI=3.700000047683716\n",
            "loading 364...0.6 -> 0.53, no AKI, min CI=1.2000000476837158, max CI=3.200000047683716\n",
            "loading 369...0.61 -> 0.6, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 4472...1.09 -> 1.13, no AKI, min CI=2.0999999046325684, max CI=3.700000047683716\n",
            "loading 2425...0.75 -> 0.75, no AKI, min CI=2.200000047683716, max CI=3.200000047683716\n",
            "loading 4475...11.03 -> 9.87, no AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 380...0.74 -> 0.63, no AKI, min CI=1.899999976158142, max CI=3.700000047683716\n",
            "loading 4476...no postop cr\n",
            "loading 382...no preop cr\n",
            "loading 4478...1.03 -> 1.23, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 383...2.32 -> 2.19, no AKI, min CI=1.7000000476837158, max CI=3.4000000953674316\n",
            "loading 385...0.64 -> 0.55, no AKI, min CI=1.899999976158142, max CI=2.799999952316284\n",
            "loading 2434...0.74 -> 0.67, no AKI, min CI=2.200000047683716, max CI=3.799999952316284\n",
            "loading 4480...1.03 -> 1.18, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 388...no preop cr\n",
            "loading 2436...0.59 -> 0.71, no AKI, min CI=2.0, max CI=3.700000047683716\n",
            "loading 4485...1.44 -> 1.27, no AKI, min CI=0.699999988079071, max CI=3.200000047683716\n",
            "loading 2433...no preop cr\n",
            "loading 4481...0.72 -> 0.84, no AKI, min CI=1.2000000476837158, max CI=3.9000000953674316\n",
            "loading 390...0.68 -> 0.67, no AKI, min CI=1.0, max CI=3.799999952316284\n",
            "loading 2443...0.93 -> 0.65, no AKI, min CI=0.8999999761581421, max CI=2.5999999046325684\n",
            "loading 397...1.9 -> 1.51, no AKI, min CI=2.4000000953674316, max CI=3.9000000953674316\n",
            "loading 4498...0.64 -> 0.57, no AKI, min CI=2.9000000953674316, max CI=3.9000000953674316\n",
            "loading 405...0.73 -> 0.59, no AKI, min CI=2.5, max CI=3.9000000953674316\n",
            "loading 2453...0.66 -> 0.72, no AKI, min CI=3.5999999046325684, max CI=3.9000000953674316\n",
            "loading 409...no postop cr\n",
            "loading 2458...0.77 -> 0.88, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 2460...0.6 -> 0.53, no AKI, min CI=1.5, max CI=3.0\n",
            "loading 4513...no preop cr\n",
            "loading 418...0.85 -> 3.73, AKI, min CI=1.7999999523162842, max CI=3.799999952316284\n",
            "loading 4515...15.26 -> 13.94, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 2472...1.01 -> 0.97, no AKI, min CI=1.2000000476837158, max CI=3.5\n",
            "loading 4522...0.8 -> 0.71, no AKI, min CI=1.399999976158142, max CI=2.5\n",
            "loading 4527...0.87 -> 0.89, no AKI, min CI=1.600000023841858, max CI=3.0\n",
            "loading 2480...0.46 -> 0.55, no AKI, min CI=2.0999999046325684, max CI=2.9000000953674316\n",
            "loading 440...0.65 -> 0.74, no AKI, min CI=1.899999976158142, max CI=3.700000047683716\n",
            "loading 441...no preop cr\n",
            "loading 4538...no preop cr\n",
            "loading 442...0.66 -> 0.84, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 4540...0.96 -> 0.96, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 2493...0.94 -> 1.88, AKI, min CI=1.399999976158142, max CI=3.4000000953674316\n",
            "loading 2495...no preop cr\n",
            "loading 448...8.28 -> 7.52, no AKI, min CI=3.0999999046325684, max CI=3.9000000953674316\n",
            "loading 2497...0.87 -> 0.81, no AKI, min CI=2.5999999046325684, max CI=3.9000000953674316\n",
            "loading 449...0.7 -> 0.73, no AKI, min CI=1.2999999523162842, max CI=3.9000000953674316\n",
            "loading 4547...no preop cr\n",
            "loading 452...0.81 -> 0.84, no AKI, min CI=1.600000023841858, max CI=2.5999999046325684\n",
            "loading 4549...0.89 -> 0.87, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 455...0.88 -> 0.87, no AKI, min CI=2.5, max CI=3.9000000953674316\n",
            "loading 2504...0.72 -> 0.61, no AKI, min CI=1.600000023841858, max CI=3.299999952316284\n",
            "loading 2507...0.73 -> 0.66, no AKI, min CI=1.2000000476837158, max CI=3.700000047683716\n",
            "loading 2508...0.84 -> 0.95, no AKI, min CI=1.0, max CI=3.5\n",
            "loading 4558...6.7 -> 3.18, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 2511...0.94 -> 0.96, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 4559...0.67 -> 0.85, no AKI, min CI=1.600000023841858, max CI=2.9000000953674316\n",
            "loading 472...0.68 -> 0.57, no AKI, min CI=2.200000047683716, max CI=3.9000000953674316\n",
            "loading 4569...no preop cr\n",
            "loading 474...0.96 -> 0.95, no AKI, min CI=1.399999976158142, max CI=3.4000000953674316\n",
            "loading 4573...0.98 -> 0.89, no AKI, min CI=1.2999999523162842, max CI=3.700000047683716\n",
            "loading 4576...8.36 -> 7.63, no AKI, min CI=1.600000023841858, max CI=3.9000000953674316\n",
            "loading 4578...1.56 -> 1.97, no AKI, min CI=3.299999952316284, max CI=3.9000000953674316\n",
            "loading 485...0.8 -> 0.73, no AKI, min CI=0.699999988079071, max CI=3.299999952316284\n",
            "loading 486...1.24 -> 1.08, no AKI, min CI=0.10000000149011612, max CI=3.9000000953674316\n",
            "loading 488...0.63 -> 0.69, no AKI, min CI=2.0999999046325684, max CI=3.799999952316284\n",
            "loading 492...0.98 -> 0.91, no AKI, min CI=1.0, max CI=3.0\n",
            "loading 2542...0.6 -> 0.65, no AKI, min CI=1.899999976158142, max CI=3.299999952316284\n",
            "loading 4597...1.24 -> 1.32, no AKI, min CI=1.5, max CI=3.299999952316284\n",
            "loading 4599...6.88 -> 5.43, no AKI, min CI=1.399999976158142, max CI=3.9000000953674316\n",
            "loading 505...1.22 -> 1.21, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 507...no ci\n",
            "loading 4604...0.54 -> 0.47, no AKI, min CI=2.200000047683716, max CI=3.9000000953674316\n",
            "loading 512...no preop cr\n",
            "loading 4609...0.71 -> 0.79, no AKI, min CI=1.600000023841858, max CI=3.9000000953674316\n",
            "loading 516...0.6 -> 0.48, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 4612...0.55 -> 0.43, no AKI, min CI=0.699999988079071, max CI=2.799999952316284\n",
            "loading 2566...1.05 -> 1.11, no AKI, min CI=2.5999999046325684, max CI=3.9000000953674316\n",
            "loading 4617...1.03 -> 0.85, no AKI, min CI=0.20000000298023224, max CI=3.4000000953674316\n",
            "loading 4620...0.98 -> 0.98, no AKI, min CI=1.2000000476837158, max CI=3.700000047683716\n",
            "loading 4621...8.56 -> 7.43, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 526...0.86 -> 1.02, no AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 530...no preop cr\n",
            "loading 2578...0.67 -> 0.71, no AKI, min CI=1.0, max CI=3.4000000953674316\n",
            "loading 533...0.94 -> 0.78, no AKI, min CI=1.0, max CI=3.9000000953674316\n",
            "loading 2583...0.71 -> 0.59, no AKI, min CI=1.7000000476837158, max CI=3.5999999046325684\n",
            "loading 4632...no preop cr\n",
            "loading 2585...no preop cr\n",
            "loading 2589...1.25 -> 0.93, no AKI, min CI=0.800000011920929, max CI=3.9000000953674316\n",
            "loading 544...no preop cr\n",
            "loading 545...6.44 -> 7.46, no AKI, min CI=3.299999952316284, max CI=3.9000000953674316\n",
            "loading 2593...12.24 -> 12.31, no AKI, min CI=1.0, max CI=3.9000000953674316\n",
            "loading 2594...1.9 -> 2.19, no AKI, min CI=1.899999976158142, max CI=3.4000000953674316\n",
            "loading 4644...no preop cr\n",
            "loading 2597...0.53 -> 0.52, no AKI, min CI=1.5, max CI=3.4000000953674316\n",
            "loading 553...0.73 -> 0.83, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 2602...0.78 -> 0.77, no AKI, min CI=1.2000000476837158, max CI=3.700000047683716\n",
            "loading 4650...19.99 -> 17.48, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 4652...0.82 -> 1.08, no AKI, min CI=1.399999976158142, max CI=3.9000000953674316\n",
            "loading 2609...no preop cr\n",
            "loading 4658...0.58 -> 0.5, no AKI, min CI=0.10000000149011612, max CI=3.9000000953674316\n",
            "loading 4660...no postop cr\n",
            "loading 564...no preop cr\n",
            "loading 566...1.09 -> 1.12, no AKI, min CI=1.600000023841858, max CI=3.9000000953674316\n",
            "loading 4662...1.04 -> 1.05, no AKI, min CI=3.700000047683716, max CI=3.9000000953674316\n",
            "loading 2613...0.75 -> 0.59, no AKI, min CI=3.0999999046325684, max CI=3.9000000953674316\n",
            "loading 4665...no preop cr\n",
            "loading 2618...0.61 -> 0.77, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 4666...0.82 -> 0.74, no AKI, min CI=2.200000047683716, max CI=3.9000000953674316\n",
            "loading 4673...1.1 -> 0.94, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 2628...0.92 -> 0.96, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 582...0.95 -> 0.85, no AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 2631...0.69 -> 0.65, no AKI, min CI=1.899999976158142, max CI=3.799999952316284\n",
            "loading 584...1.19 -> 1.05, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 2635...no preop cr\n",
            "loading 4684...0.68 -> 0.73, no AKI, min CI=0.10000000149011612, max CI=3.9000000953674316\n",
            "loading 2637...no preop cr\n",
            "loading 2641...1.11 -> 1.07, no AKI, min CI=1.7999999523162842, max CI=3.299999952316284\n",
            "loading 599...0.9 -> 1.0, no AKI, min CI=1.5, max CI=2.700000047683716\n",
            "loading 4695...no postop cr\n",
            "loading 2648...0.8 -> 0.79, no AKI, min CI=1.2000000476837158, max CI=2.799999952316284\n",
            "loading 2654...1.22 -> 1.13, no AKI, min CI=1.100000023841858, max CI=3.9000000953674316\n",
            "loading 4703...no postop cr\n",
            "loading 610...0.8 -> 0.79, no AKI, min CI=2.5, max CI=3.0999999046325684\n",
            "loading 2662...no preop cr\n",
            "loading 2667...0.59 -> 0.57, no AKI, min CI=0.10000000149011612, max CI=3.9000000953674316\n",
            "loading 4718...0.51 -> 0.66, no AKI, min CI=1.5, max CI=3.200000047683716\n",
            "loading 624...0.76 -> 0.81, no AKI, min CI=1.2000000476837158, max CI=3.9000000953674316\n",
            "loading 4721...0.71 -> 0.65, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 629...0.82 -> 0.98, no AKI, min CI=0.4000000059604645, max CI=3.9000000953674316\n",
            "loading 631...8.24 -> 6.95, no AKI, min CI=2.700000047683716, max CI=3.9000000953674316\n",
            "loading 634...0.84 -> 0.88, no AKI, min CI=2.200000047683716, max CI=3.700000047683716\n",
            "loading 638...no postop cr\n",
            "loading 641...0.82 -> 0.91, no AKI, min CI=1.2000000476837158, max CI=3.9000000953674316\n",
            "loading 2690...0.87 -> 0.9, no AKI, min CI=1.100000023841858, max CI=3.5999999046325684\n",
            "loading 4743...no preop cr\n",
            "loading 4744...5.85 -> 5.61, no AKI, min CI=2.5, max CI=3.9000000953674316\n",
            "loading 4746...0.89 -> 0.99, no AKI, min CI=0.20000000298023224, max CI=3.9000000953674316\n",
            "loading 2699...no preop cr\n",
            "loading 2700...3.85 -> 3.06, no AKI, min CI=3.0, max CI=3.9000000953674316\n",
            "loading 654...0.72 -> 0.6, no AKI, min CI=2.799999952316284, max CI=3.9000000953674316\n",
            "loading 655...0.82 -> 0.75, no AKI, min CI=1.7000000476837158, max CI=3.700000047683716\n",
            "loading 2716...5.67 -> 4.36, no AKI, min CI=2.299999952316284, max CI=3.700000047683716\n",
            "loading 2720...0.81 -> 0.81, no AKI, min CI=3.0999999046325684, max CI=3.9000000953674316\n",
            "loading 2724...0.52 -> 0.41, no AKI, min CI=1.2000000476837158, max CI=2.799999952316284\n",
            "loading 4778...0.76 -> 0.82, no AKI, min CI=0.8999999761581421, max CI=3.9000000953674316\n",
            "loading 2732...0.73 -> 0.61, no AKI, min CI=1.7999999523162842, max CI=3.700000047683716\n",
            "loading 4781...1.3 -> 1.05, no AKI, min CI=1.600000023841858, max CI=3.799999952316284\n",
            "loading 4780...0.84 -> 0.84, no AKI, min CI=1.5, max CI=2.5\n",
            "loading 2733...0.61 -> 0.57, no AKI, min CI=0.30000001192092896, max CI=2.799999952316284\n",
            "loading 2735...1.17 -> 1.1, no AKI, min CI=1.399999976158142, max CI=3.200000047683716\n",
            "loading 2738...0.73 -> 1.27, AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 691...11.65 -> 10.3, no AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 4789...0.72 -> 0.68, no AKI, min CI=2.0, max CI=3.5\n",
            "loading 4793...no preop cr\n",
            "loading 2746...0.7 -> 0.51, no AKI, min CI=1.7000000476837158, max CI=2.5\n",
            "loading 4798...0.67 -> 0.56, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 2751...1.01 -> 1.02, no AKI, min CI=1.2999999523162842, max CI=3.9000000953674316\n",
            "loading 4800...no ci\n",
            "loading 706...no postop cr\n",
            "loading 711...1.12 -> 1.11, no AKI, min CI=1.100000023841858, max CI=3.9000000953674316\n",
            "loading 4808...no preop cr\n",
            "loading 4809...0.94 -> 0.9, no AKI, min CI=1.600000023841858, max CI=3.9000000953674316\n",
            "loading 2761...no postop cr\n",
            "loading 2765...no preop cr\n",
            "loading 2766...no preop cr\n",
            "loading 2769...no preop cr\n",
            "loading 722...no preop cr\n",
            "loading 4820...0.89 -> 1.08, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 725...no preop cr\n",
            "loading 2778...1.2 -> 2.27, AKI, min CI=1.2999999523162842, max CI=3.799999952316284\n",
            "loading 2779...0.87 -> 0.81, no AKI, min CI=2.4000000953674316, max CI=3.5999999046325684\n",
            "loading 2781...0.85 -> 0.82, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 734...1.19 -> 1.09, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 4834...0.79 -> 0.79, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 2788...0.87 -> 1.02, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 741...0.86 -> 0.76, no AKI, min CI=1.600000023841858, max CI=3.9000000953674316\n",
            "loading 4837...1.09 -> 0.83, no AKI, min CI=1.7999999523162842, max CI=3.0\n",
            "loading 740...0.9 -> 0.95, no AKI, min CI=2.4000000953674316, max CI=2.5999999046325684\n",
            "loading 742...1.06 -> 1.03, no AKI, min CI=1.2000000476837158, max CI=3.9000000953674316\n",
            "loading 749...no postop cr\n",
            "loading 751...1.03 -> 1.06, no AKI, min CI=2.200000047683716, max CI=3.9000000953674316\n",
            "loading 4851...0.44 -> 0.47, no AKI, min CI=0.8999999761581421, max CI=3.9000000953674316\n",
            "loading 4857...0.92 -> 0.78, no AKI, min CI=1.2000000476837158, max CI=3.9000000953674316\n",
            "loading 4858...0.61 -> 0.53, no AKI, min CI=2.200000047683716, max CI=3.9000000953674316\n",
            "loading 768...0.62 -> 0.66, no AKI, min CI=1.899999976158142, max CI=2.5\n",
            "loading 772...0.76 -> 0.9, no AKI, min CI=2.0, max CI=3.799999952316284\n",
            "loading 773...0.63 -> 0.71, no AKI, min CI=1.600000023841858, max CI=3.0999999046325684\n",
            "loading 775...0.59 -> 1.04, AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 4871...1.16 -> 1.07, no AKI, min CI=0.8999999761581421, max CI=2.9000000953674316\n",
            "loading 777...0.57 -> 0.66, no AKI, min CI=2.200000047683716, max CI=3.9000000953674316\n",
            "loading 779...0.6 -> 0.52, no AKI, min CI=3.200000047683716, max CI=3.9000000953674316\n",
            "loading 781...no postop cr\n",
            "loading 4879...1.05 -> 1.03, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 783...1.51 -> 1.58, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 4883...1.01 -> 0.95, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 787...no preop cr\n",
            "loading 2838...0.8 -> 0.92, no AKI, min CI=0.10000000149011612, max CI=3.9000000953674316\n",
            "loading 2839...0.95 -> 0.88, no AKI, min CI=2.0, max CI=3.0\n",
            "loading 792...no preop cr\n",
            "loading 794...0.67 -> 0.8, no AKI, min CI=2.799999952316284, max CI=3.9000000953674316\n",
            "loading 2845...0.95 -> 0.87, no AKI, min CI=1.899999976158142, max CI=3.799999952316284\n",
            "loading 2849...no preop cr\n",
            "loading 2850...0.93 -> 1.07, no AKI, min CI=1.7000000476837158, max CI=3.700000047683716\n",
            "loading 4897...0.7 -> 1.18, AKI, min CI=1.899999976158142, max CI=2.9000000953674316\n",
            "loading 4903...no ci\n",
            "loading 4904...0.7 -> 0.71, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 2858...1.05 -> 1.08, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 2860...5.94 -> 6.51, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 814...0.89 -> 0.77, no AKI, min CI=0.10000000149011612, max CI=3.4000000953674316\n",
            "loading 4913...1.29 -> 0.75, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 4914...0.65 -> 0.65, no AKI, min CI=1.100000023841858, max CI=3.5\n",
            "loading 819...0.69 -> 0.71, no AKI, min CI=1.2999999523162842, max CI=2.700000047683716\n",
            "loading 818...0.92 -> 1.7, AKI, min CI=1.7999999523162842, max CI=3.5\n",
            "loading 822...0.57 -> 0.52, no AKI, min CI=2.200000047683716, max CI=3.9000000953674316\n",
            "loading 2872...0.91 -> 0.9, no AKI, min CI=1.2999999523162842, max CI=3.9000000953674316\n",
            "loading 831...no postop cr\n",
            "loading 4929...0.73 -> 0.62, no AKI, min CI=2.5, max CI=3.9000000953674316\n",
            "loading 4932...0.72 -> 0.67, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 4933...0.66 -> 0.6, no AKI, min CI=0.5, max CI=3.9000000953674316\n",
            "loading 4936...0.79 -> 0.79, no AKI, min CI=1.7000000476837158, max CI=3.299999952316284\n",
            "loading 843...0.91 -> 1.03, no AKI, min CI=0.10000000149011612, max CI=3.9000000953674316\n",
            "loading 844...0.82 -> 0.69, no AKI, min CI=0.10000000149011612, max CI=3.9000000953674316\n",
            "loading 4941...no preop cr\n",
            "loading 2895...0.61 -> 0.78, no AKI, min CI=1.7999999523162842, max CI=2.700000047683716\n",
            "loading 848...0.57 -> 0.53, no AKI, min CI=0.800000011920929, max CI=3.5\n",
            "loading 4946...0.5 -> 0.31, no AKI, min CI=1.100000023841858, max CI=3.700000047683716\n",
            "loading 4951...1.13 -> 1.06, no AKI, min CI=1.0, max CI=3.9000000953674316\n",
            "loading 4953...1.0 -> 0.75, no AKI, min CI=1.5, max CI=3.5\n",
            "loading 2906...0.85 -> 0.86, no AKI, min CI=2.200000047683716, max CI=2.700000047683716\n",
            "loading 860...0.86 -> 0.88, no AKI, min CI=1.600000023841858, max CI=3.799999952316284\n",
            "loading 4957...0.86 -> 0.87, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 2911...0.71 -> 0.75, no AKI, min CI=2.9000000953674316, max CI=3.9000000953674316\n",
            "loading 4959...no ci\n",
            "loading 866...no preop cr\n",
            "loading 2914...0.95 -> 0.97, no AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 2916...0.82 -> 0.84, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 870...9.79 -> 9.32, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 4966...1.14 -> 1.07, no AKI, min CI=0.800000011920929, max CI=3.799999952316284\n",
            "loading 872...0.86 -> 0.87, no AKI, min CI=3.200000047683716, max CI=3.9000000953674316\n",
            "loading 2922...1.35 -> 1.48, no AKI, min CI=1.100000023841858, max CI=2.5999999046325684\n",
            "loading 2930...0.66 -> 0.73, no AKI, min CI=2.200000047683716, max CI=3.9000000953674316\n",
            "loading 2933...0.67 -> 0.68, no AKI, min CI=1.7999999523162842, max CI=3.200000047683716\n",
            "loading 4982...0.72 -> 0.67, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 885...0.99 -> 0.81, no AKI, min CI=3.0999999046325684, max CI=3.799999952316284\n",
            "loading 886...0.58 -> 0.5, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 2935...1.1 -> 1.28, no AKI, min CI=1.600000023841858, max CI=3.200000047683716\n",
            "loading 4987...0.99 -> 0.77, no AKI, min CI=0.699999988079071, max CI=3.200000047683716\n",
            "loading 2943...0.78 -> 0.65, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 2945...0.5 -> 0.9, AKI, min CI=1.0, max CI=3.9000000953674316\n",
            "loading 4997...no preop cr\n",
            "loading 2952...0.96 -> 1.03, no AKI, min CI=0.10000000149011612, max CI=3.9000000953674316\n",
            "loading 907...no preop cr\n",
            "loading 5004...0.87 -> 0.83, no AKI, min CI=1.399999976158142, max CI=3.5\n",
            "loading 2958...0.76 -> 0.57, no AKI, min CI=1.5, max CI=3.0999999046325684\n",
            "loading 917...1.18 -> 1.27, no AKI, min CI=1.7000000476837158, max CI=3.799999952316284\n",
            "loading 2966...0.69 -> 0.86, no AKI, min CI=2.0999999046325684, max CI=3.799999952316284\n",
            "loading 5018...1.03 -> 1.13, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 2971...0.64 -> 0.72, no AKI, min CI=3.0, max CI=3.9000000953674316\n",
            "loading 2972...8.43 -> 7.68, no AKI, min CI=0.10000000149011612, max CI=3.9000000953674316\n",
            "loading 5019...0.99 -> 0.91, no AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 2974...0.83 -> 0.82, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 931...3.02 -> 4.13, no AKI, min CI=2.5999999046325684, max CI=3.9000000953674316\n",
            "loading 933...0.72 -> 0.85, no AKI, min CI=1.5, max CI=3.799999952316284\n",
            "loading 5029...0.64 -> 0.65, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 936...no preop cr\n",
            "loading 5033...no postop cr\n",
            "loading 937...0.66 -> 0.55, no AKI, min CI=1.0, max CI=3.9000000953674316\n",
            "loading 5035...2.61 -> 2.11, no AKI, min CI=1.399999976158142, max CI=3.200000047683716\n",
            "loading 945...0.61 -> 0.65, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 946...0.64 -> 0.55, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 5045...0.75 -> 0.65, no AKI, min CI=2.5999999046325684, max CI=3.9000000953674316\n",
            "loading 2998...no preop cr\n",
            "loading 3000...1.21 -> 1.32, no AKI, min CI=1.2999999523162842, max CI=3.799999952316284\n",
            "loading 958...0.81 -> 0.76, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 5059...0.74 -> 1.34, AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 5060...0.68 -> 1.07, AKI, min CI=0.8999999761581421, max CI=3.9000000953674316\n",
            "loading 3015...1.06 -> 0.85, no AKI, min CI=2.0999999046325684, max CI=3.5999999046325684\n",
            "loading 969...no preop cr\n",
            "loading 972...1.34 -> 1.64, no AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 3020...0.74 -> 0.58, no AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 976...0.79 -> 0.88, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 3024...0.68 -> 0.52, no AKI, min CI=1.399999976158142, max CI=3.9000000953674316\n",
            "loading 980...0.92 -> 1.01, no AKI, min CI=1.399999976158142, max CI=3.9000000953674316\n",
            "loading 3028...no preop cr\n",
            "loading 5079...0.73 -> 0.59, no AKI, min CI=0.4000000059604645, max CI=3.299999952316284\n",
            "loading 985...no postop cr\n",
            "loading 3034...0.59 -> 0.53, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 986...6.51 -> 5.82, no AKI, min CI=2.700000047683716, max CI=3.9000000953674316\n",
            "loading 5084...0.46 -> 0.47, no AKI, min CI=1.399999976158142, max CI=2.9000000953674316\n",
            "loading 3037...0.82 -> 0.7, no AKI, min CI=1.0, max CI=3.9000000953674316\n",
            "loading 3040...0.68 -> 0.63, no AKI, min CI=1.2000000476837158, max CI=3.0999999046325684\n",
            "loading 5090...no preop cr\n",
            "loading 5091...6.46 -> 6.45, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 5092...0.68 -> 0.65, no AKI, min CI=2.4000000953674316, max CI=3.9000000953674316\n",
            "loading 5093...0.64 -> 0.67, no AKI, min CI=0.800000011920929, max CI=2.799999952316284\n",
            "loading 3047...0.76 -> 0.77, no AKI, min CI=0.8999999761581421, max CI=3.799999952316284\n",
            "loading 3050...0.95 -> 0.78, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 5099...0.99 -> 0.86, no AKI, min CI=1.100000023841858, max CI=3.799999952316284\n",
            "loading 1004...0.67 -> 0.57, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 3053...no postop cr\n",
            "loading 5106...no preop cr\n",
            "loading 5107...0.58 -> 0.82, no AKI, min CI=0.8999999761581421, max CI=3.9000000953674316\n",
            "loading 3058...0.97 -> 0.79, no AKI, min CI=0.4000000059604645, max CI=3.799999952316284\n",
            "loading 1013...no preop cr\n",
            "loading 5109...0.59 -> 0.6, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 5111...0.88 -> 0.85, no AKI, min CI=1.600000023841858, max CI=2.9000000953674316\n",
            "loading 5112...no postop cr\n",
            "loading 1018...0.85 -> 0.91, no AKI, min CI=0.10000000149011612, max CI=3.9000000953674316\n",
            "loading 1025...0.53 -> 0.39, no AKI, min CI=1.399999976158142, max CI=3.200000047683716\n",
            "loading 1027...0.85 -> 0.64, no AKI, min CI=1.7999999523162842, max CI=3.5\n",
            "loading 1029...8.18 -> 1.57, no AKI, min CI=3.299999952316284, max CI=3.9000000953674316\n",
            "loading 1030...0.87 -> 0.67, no AKI, min CI=1.899999976158142, max CI=2.5999999046325684\n",
            "loading 5127...3.24 -> 2.82, no AKI, min CI=1.100000023841858, max CI=3.9000000953674316\n",
            "loading 1032...0.47 -> 0.47, no AKI, min CI=1.600000023841858, max CI=3.9000000953674316\n",
            "loading 3080...no preop cr\n",
            "loading 1034...0.66 -> 0.67, no AKI, min CI=1.600000023841858, max CI=3.9000000953674316\n",
            "loading 1035...1.04 -> 1.04, no AKI, min CI=1.7999999523162842, max CI=3.700000047683716\n",
            "loading 3085...0.94 -> 0.9, no AKI, min CI=0.800000011920929, max CI=2.9000000953674316\n",
            "loading 1038...no preop cr\n",
            "loading 5135...0.83 -> 0.75, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 3088...0.86 -> 0.86, no AKI, min CI=1.2999999523162842, max CI=3.9000000953674316\n",
            "loading 1041...0.67 -> 0.58, no AKI, min CI=2.299999952316284, max CI=3.700000047683716\n",
            "loading 3090...0.77 -> 1.02, no AKI, min CI=1.600000023841858, max CI=2.299999952316284\n",
            "loading 1044...0.96 -> 0.88, no AKI, min CI=1.399999976158142, max CI=3.0\n",
            "loading 5141...no preop cr\n",
            "loading 3093...no postop cr\n",
            "loading 3094...0.93 -> 0.98, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 3096...0.89 -> 0.94, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 5140...no preop cr\n",
            "loading 5142...0.96 -> 1.11, no AKI, min CI=2.200000047683716, max CI=3.9000000953674316\n",
            "loading 3098...no preop cr\n",
            "loading 1051...1.05 -> 0.9, no AKI, min CI=1.2999999523162842, max CI=2.5999999046325684\n",
            "loading 3103...1.01 -> 1.15, no AKI, min CI=2.9000000953674316, max CI=3.9000000953674316\n",
            "loading 1056...8.4 -> 7.36, no AKI, min CI=2.4000000953674316, max CI=3.9000000953674316\n",
            "loading 5153...0.92 -> 0.94, no AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 3106...1.14 -> 0.93, no AKI, min CI=2.799999952316284, max CI=3.700000047683716\n",
            "loading 3107...0.59 -> 0.54, no AKI, min CI=2.299999952316284, max CI=3.799999952316284\n",
            "loading 1061...15.15 -> 13.71, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 3113...0.92 -> 1.0, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 5162...0.91 -> 0.93, no AKI, min CI=0.699999988079071, max CI=2.700000047683716\n",
            "loading 3122...0.74 -> 0.68, no AKI, min CI=1.5, max CI=3.299999952316284\n",
            "loading 1076...0.79 -> 0.91, no AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 3125...0.74 -> 0.72, no AKI, min CI=1.0, max CI=2.4000000953674316\n",
            "loading 3126...1.15 -> 1.2, no AKI, min CI=2.0999999046325684, max CI=3.700000047683716\n",
            "loading 1078...0.95 -> 1.15, no AKI, min CI=0.800000011920929, max CI=3.799999952316284\n",
            "loading 5175...no postop cr\n",
            "loading 3130...no preop cr\n",
            "loading 1083...no preop cr\n",
            "loading 5182...no postop cr\n",
            "loading 1095...0.95 -> 0.88, no AKI, min CI=2.4000000953674316, max CI=3.9000000953674316\n",
            "loading 1097...0.79 -> 0.72, no AKI, min CI=1.899999976158142, max CI=3.0999999046325684\n",
            "loading 3145...no postop cr\n",
            "loading 3148...0.97 -> 1.05, no AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 3149...0.98 -> 0.98, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 3150...0.82 -> 0.69, no AKI, min CI=2.5999999046325684, max CI=3.700000047683716\n",
            "loading 5200...14.37 -> 11.1, no AKI, min CI=1.600000023841858, max CI=3.9000000953674316\n",
            "loading 3159...0.84 -> 0.85, no AKI, min CI=1.2999999523162842, max CI=3.9000000953674316\n",
            "loading 1113...0.44 -> 0.56, no AKI, min CI=1.399999976158142, max CI=3.200000047683716\n",
            "loading 3170...0.68 -> 0.52, no AKI, min CI=1.399999976158142, max CI=3.799999952316284\n",
            "loading 1123...0.73 -> 0.79, no AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 3172...0.84 -> 0.86, no AKI, min CI=2.5, max CI=3.9000000953674316\n",
            "loading 1124...1.03 -> 1.14, no AKI, min CI=1.0, max CI=3.799999952316284\n",
            "loading 5222...5.46 -> 3.49, no AKI, min CI=0.4000000059604645, max CI=3.700000047683716\n",
            "loading 3176...1.26 -> 1.33, no AKI, min CI=2.0, max CI=3.799999952316284\n",
            "loading 5226...1.56 -> 1.34, no AKI, min CI=2.5999999046325684, max CI=3.9000000953674316\n",
            "loading 1132...no postop cr\n",
            "loading 1133...0.77 -> 0.85, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 3188...no postop cr\n",
            "loading 3189...no preop cr\n",
            "loading 3193...no postop cr\n",
            "loading 3196...0.55 -> 0.52, no AKI, min CI=1.7999999523162842, max CI=3.799999952316284\n",
            "loading 5251...13.12 -> 13.43, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 3207...no postop cr\n",
            "loading 5258...0.81 -> 0.66, no AKI, min CI=2.0, max CI=3.5\n",
            "loading 5259...0.31 -> 0.38, no AKI, min CI=2.9000000953674316, max CI=3.9000000953674316\n",
            "loading 1165...0.99 -> 1.02, no AKI, min CI=0.8999999761581421, max CI=3.9000000953674316\n",
            "loading 1166...8.62 -> 1.7, no AKI, min CI=2.700000047683716, max CI=3.9000000953674316\n",
            "loading 1172...0.52 -> 0.54, no AKI, min CI=1.100000023841858, max CI=2.5999999046325684\n",
            "loading 5269...1.43 -> 1.78, no AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 3223...0.89 -> 0.97, no AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 5278...0.85 -> 0.78, no AKI, min CI=1.2000000476837158, max CI=3.0999999046325684\n",
            "loading 3232...no preop cr\n",
            "loading 1185...0.79 -> 0.82, no AKI, min CI=1.899999976158142, max CI=3.0999999046325684\n",
            "loading 5282...1.38 -> 1.36, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 5283...0.58 -> 0.59, no AKI, min CI=2.9000000953674316, max CI=3.9000000953674316\n",
            "loading 3235...no preop cr\n",
            "loading 1191...0.62 -> 1.06, AKI, min CI=0.800000011920929, max CI=3.9000000953674316\n",
            "loading 1193...0.69 -> 0.68, no AKI, min CI=1.7999999523162842, max CI=3.5\n",
            "loading 1194...0.6 -> 0.6, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 5290...no ci\n",
            "loading 5295...0.53 -> 0.47, no AKI, min CI=2.200000047683716, max CI=3.9000000953674316\n",
            "loading 3248...12.03 -> 10.59, no AKI, min CI=1.2999999523162842, max CI=3.9000000953674316\n",
            "loading 5296...0.86 -> 0.75, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 5299...no preop cr\n",
            "loading 3255...0.91 -> 0.72, no AKI, min CI=1.600000023841858, max CI=3.9000000953674316\n",
            "loading 1208...0.83 -> 0.87, no AKI, min CI=1.2999999523162842, max CI=3.9000000953674316\n",
            "loading 5304...0.76 -> 1.06, no AKI, min CI=1.0, max CI=3.9000000953674316\n",
            "loading 3258...0.31 -> 0.37, no AKI, min CI=1.100000023841858, max CI=3.700000047683716\n",
            "loading 1213...1.18 -> 1.12, no AKI, min CI=1.399999976158142, max CI=3.5999999046325684\n",
            "loading 5311...0.52 -> 0.53, no AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 1215...0.66 -> 0.58, no AKI, min CI=1.399999976158142, max CI=3.9000000953674316\n",
            "loading 1219...0.76 -> 0.73, no AKI, min CI=1.600000023841858, max CI=2.4000000953674316\n",
            "loading 3267...0.75 -> 0.6, no AKI, min CI=1.7000000476837158, max CI=3.799999952316284\n",
            "loading 1221...0.85 -> 1.0, no AKI, min CI=0.8999999761581421, max CI=2.5\n",
            "loading 5317...0.93 -> 1.08, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 3270...0.66 -> 0.69, no AKI, min CI=1.399999976158142, max CI=3.9000000953674316\n",
            "loading 5319...1.02 -> 0.75, no AKI, min CI=1.2000000476837158, max CI=3.9000000953674316\n",
            "loading 1225...0.74 -> 0.57, no AKI, min CI=1.100000023841858, max CI=3.5999999046325684\n",
            "loading 3275...0.69 -> 0.66, no AKI, min CI=1.600000023841858, max CI=3.9000000953674316\n",
            "loading 5323...0.83 -> 0.73, no AKI, min CI=1.7999999523162842, max CI=3.4000000953674316\n",
            "loading 3279...0.79 -> 0.83, no AKI, min CI=1.100000023841858, max CI=3.9000000953674316\n",
            "loading 5327...0.74 -> 0.65, no AKI, min CI=2.4000000953674316, max CI=3.799999952316284\n",
            "loading 1233...0.84 -> 0.81, no AKI, min CI=1.7000000476837158, max CI=2.799999952316284\n",
            "loading 5328...0.65 -> 0.95, no AKI, min CI=1.2000000476837158, max CI=3.9000000953674316\n",
            "loading 1237...0.88 -> 0.87, no AKI, min CI=2.700000047683716, max CI=3.9000000953674316\n",
            "loading 1239...1.0 -> 0.85, no AKI, min CI=1.2999999523162842, max CI=3.9000000953674316\n",
            "loading 3287...0.73 -> 0.79, no AKI, min CI=2.5999999046325684, max CI=3.9000000953674316\n",
            "loading 3291...0.54 -> 0.69, no AKI, min CI=0.4000000059604645, max CI=3.9000000953674316\n",
            "loading 3296...0.63 -> 0.58, no AKI, min CI=1.5, max CI=3.299999952316284\n",
            "loading 1249...no preop cr\n",
            "loading 5346...1.26 -> 1.4, no AKI, min CI=1.600000023841858, max CI=3.9000000953674316\n",
            "loading 3299...0.91 -> 0.61, no AKI, min CI=1.399999976158142, max CI=3.799999952316284\n",
            "loading 3300...5.28 -> 4.61, no AKI, min CI=1.399999976158142, max CI=3.9000000953674316\n",
            "loading 5349...8.08 -> 5.5, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 5350...no ci\n",
            "loading 1256...0.95 -> 1.07, no AKI, min CI=1.100000023841858, max CI=3.4000000953674316\n",
            "loading 3310...no preop cr\n",
            "loading 3311...0.78 -> 1.74, AKI, min CI=0.10000000149011612, max CI=3.9000000953674316\n",
            "loading 5360...1.37 -> 1.04, no AKI, min CI=0.5, max CI=3.5999999046325684\n",
            "loading 1267...0.93 -> 0.96, no AKI, min CI=1.2000000476837158, max CI=3.200000047683716\n",
            "loading 3317...0.65 -> 0.71, no AKI, min CI=2.200000047683716, max CI=3.0999999046325684\n",
            "loading 3325...0.76 -> 0.63, no AKI, min CI=2.5, max CI=3.9000000953674316\n",
            "loading 3328...0.58 -> 1.84, AKI, min CI=1.2000000476837158, max CI=3.9000000953674316\n",
            "loading 1281...0.5 -> 0.42, no AKI, min CI=2.5, max CI=3.9000000953674316\n",
            "loading 3329...0.84 -> 0.98, no AKI, min CI=1.2999999523162842, max CI=3.200000047683716\n",
            "loading 3330...0.94 -> 0.89, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 1286...0.66 -> 0.44, no AKI, min CI=0.20000000298023224, max CI=3.5\n",
            "loading 1290...no preop cr\n",
            "loading 1292...0.81 -> 1.02, no AKI, min CI=0.800000011920929, max CI=3.9000000953674316\n",
            "loading 5396...1.06 -> 3.95, AKI, min CI=0.800000011920929, max CI=3.9000000953674316\n",
            "loading 1301...0.85 -> 0.9, no AKI, min CI=2.0, max CI=3.799999952316284\n",
            "loading 3348...0.58 -> 0.51, no AKI, min CI=1.5, max CI=2.5\n",
            "loading 1303...0.66 -> 0.69, no AKI, min CI=1.7000000476837158, max CI=2.5\n",
            "loading 1302...0.65 -> 0.58, no AKI, min CI=1.7000000476837158, max CI=3.700000047683716\n",
            "loading 1307...0.89 -> 0.99, no AKI, min CI=1.399999976158142, max CI=3.799999952316284\n",
            "loading 5403...0.68 -> 0.91, no AKI, min CI=0.8999999761581421, max CI=3.799999952316284\n",
            "loading 5404...0.59 -> 0.66, no AKI, min CI=1.100000023841858, max CI=3.0999999046325684\n",
            "loading 3363...no preop cr\n",
            "loading 5411...0.75 -> 0.73, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 3364...0.68 -> 0.6, no AKI, min CI=2.299999952316284, max CI=3.5999999046325684\n",
            "loading 1321...1.48 -> 1.56, no AKI, min CI=2.0999999046325684, max CI=3.5999999046325684\n",
            "loading 1323...no postop cr\n",
            "loading 1325...0.9 -> 1.0, no AKI, min CI=1.2999999523162842, max CI=3.9000000953674316\n",
            "loading 5422...6.56 -> 4.86, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 1327...1.19 -> 1.63, no AKI, min CI=0.20000000298023224, max CI=3.9000000953674316\n",
            "loading 5425...6.9 -> 6.55, no AKI, min CI=2.4000000953674316, max CI=3.9000000953674316\n",
            "loading 5428...1.13 -> 0.87, no AKI, min CI=0.800000011920929, max CI=3.9000000953674316\n",
            "loading 1335...0.74 -> 0.65, no AKI, min CI=1.100000023841858, max CI=2.5\n",
            "loading 3385...1.03 -> 1.13, no AKI, min CI=3.0, max CI=3.9000000953674316\n",
            "loading 1339...no postop cr\n",
            "loading 3391...1.04 -> 0.97, no AKI, min CI=2.5999999046325684, max CI=3.9000000953674316\n",
            "loading 5439...0.85 -> 0.78, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 1346...1.04 -> 0.88, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 3394...0.74 -> 0.71, no AKI, min CI=1.600000023841858, max CI=3.9000000953674316\n",
            "loading 5442...no postop cr\n",
            "loading 1350...0.8 -> 0.73, no AKI, min CI=2.9000000953674316, max CI=3.9000000953674316\n",
            "loading 5453...1.26 -> 1.17, no AKI, min CI=1.7000000476837158, max CI=3.4000000953674316\n",
            "loading 3407...no postop cr\n",
            "loading 5460...0.71 -> 0.64, no AKI, min CI=1.600000023841858, max CI=3.299999952316284\n",
            "loading 3413...1.66 -> 1.83, no AKI, min CI=3.200000047683716, max CI=3.9000000953674316\n",
            "loading 5463...0.56 -> 0.54, no AKI, min CI=1.7999999523162842, max CI=3.700000047683716\n",
            "loading 5464...0.6 -> 0.53, no AKI, min CI=1.7999999523162842, max CI=2.700000047683716\n",
            "loading 3418...0.81 -> 0.94, no AKI, min CI=2.5, max CI=3.9000000953674316\n",
            "loading 5467...0.65 -> 0.65, no AKI, min CI=1.0, max CI=3.9000000953674316\n",
            "loading 1376...0.64 -> 0.51, no AKI, min CI=1.899999976158142, max CI=3.4000000953674316\n",
            "loading 5472...4.02 -> 3.51, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 3426...0.48 -> 0.56, no AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 5480...0.86 -> 0.68, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 3438...0.55 -> 1.5, AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 5490...0.73 -> 0.57, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 1396...0.84 -> 0.83, no AKI, min CI=1.2000000476837158, max CI=3.799999952316284\n",
            "loading 1397...2.02 -> 2.43, no AKI, min CI=0.5, max CI=3.5999999046325684\n",
            "loading 3447...0.42 -> 0.79, AKI, min CI=1.600000023841858, max CI=3.9000000953674316\n",
            "loading 3449...0.56 -> 0.43, no AKI, min CI=1.100000023841858, max CI=2.700000047683716\n",
            "loading 1407...0.87 -> 2.34, AKI, min CI=0.10000000149011612, max CI=3.9000000953674316\n",
            "loading 1410...no ci\n",
            "loading 3458...0.77 -> 0.92, no AKI, min CI=0.10000000149011612, max CI=3.799999952316284\n",
            "loading 5508...0.93 -> 0.81, no AKI, min CI=1.600000023841858, max CI=3.9000000953674316\n",
            "loading 5509...0.68 -> 0.76, no AKI, min CI=0.8999999761581421, max CI=3.9000000953674316\n",
            "loading 3462...0.63 -> 0.62, no AKI, min CI=2.9000000953674316, max CI=3.9000000953674316\n",
            "loading 5517...0.82 -> 0.72, no AKI, min CI=1.399999976158142, max CI=3.799999952316284\n",
            "loading 5520...no preop cr\n",
            "loading 3475...0.52 -> 0.46, no AKI, min CI=1.2999999523162842, max CI=3.0999999046325684\n",
            "loading 3476...0.56 -> 0.55, no AKI, min CI=1.100000023841858, max CI=3.4000000953674316\n",
            "loading 3478...8.17 -> 6.95, no AKI, min CI=0.800000011920929, max CI=3.9000000953674316\n",
            "loading 1432...0.81 -> 0.7, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 5534...no postop cr\n",
            "loading 5536...0.62 -> 0.56, no AKI, min CI=1.600000023841858, max CI=3.200000047683716\n",
            "loading 3488...1.01 -> 0.95, no AKI, min CI=1.600000023841858, max CI=3.9000000953674316\n",
            "loading 1444...0.89 -> 0.96, no AKI, min CI=2.5, max CI=3.799999952316284\n",
            "loading 5546...1.24 -> 1.06, no AKI, min CI=3.299999952316284, max CI=3.9000000953674316\n",
            "loading 3501...0.61 -> 0.74, no AKI, min CI=1.899999976158142, max CI=3.4000000953674316\n",
            "loading 1454...1.19 -> 1.19, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 3503...0.81 -> 0.8, no AKI, min CI=1.600000023841858, max CI=3.0999999046325684\n",
            "loading 5550...no ci\n",
            "loading 3505...4.95 -> 4.65, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 5552...0.73 -> 0.69, no AKI, min CI=1.7999999523162842, max CI=2.799999952316284\n",
            "loading 5556...0.81 -> 0.81, no AKI, min CI=0.699999988079071, max CI=3.5\n",
            "loading 1468...no preop cr\n",
            "loading 3516...4.68 -> 3.85, no AKI, min CI=3.5999999046325684, max CI=3.9000000953674316\n",
            "loading 1469...0.63 -> 0.5, no AKI, min CI=2.5, max CI=3.9000000953674316\n",
            "loading 1473...0.97 -> 0.87, no AKI, min CI=1.399999976158142, max CI=2.299999952316284\n",
            "loading 3521...0.75 -> 0.76, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 1475...0.9 -> 0.83, no AKI, min CI=1.2000000476837158, max CI=3.0\n",
            "loading 3524...0.55 -> 0.52, no AKI, min CI=3.4000000953674316, max CI=3.9000000953674316\n",
            "loading 1478...no preop cr\n",
            "loading 5574...0.58 -> 0.51, no AKI, min CI=0.699999988079071, max CI=2.9000000953674316\n",
            "loading 3529...0.8 -> 0.71, no AKI, min CI=1.2999999523162842, max CI=3.9000000953674316\n",
            "loading 5578...0.75 -> 0.72, no AKI, min CI=2.4000000953674316, max CI=3.9000000953674316\n",
            "loading 1482...0.57 -> 0.48, no AKI, min CI=1.7999999523162842, max CI=3.5999999046325684\n",
            "loading 3531...no preop cr\n",
            "loading 1485...no postop cr\n",
            "loading 3533...3.82 -> 2.86, no AKI, min CI=0.10000000149011612, max CI=3.9000000953674316\n",
            "loading 1488...1.23 -> 1.43, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 1489...0.74 -> 0.66, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 3537...0.41 -> 0.6, no AKI, min CI=2.200000047683716, max CI=3.9000000953674316\n",
            "loading 5585...0.71 -> 0.64, no AKI, min CI=1.0, max CI=2.5999999046325684\n",
            "loading 3545...0.64 -> 0.54, no AKI, min CI=0.8999999761581421, max CI=3.9000000953674316\n",
            "loading 3546...0.71 -> 0.57, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 5595...1.16 -> 1.02, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 3550...0.42 -> 0.59, no AKI, min CI=1.2999999523162842, max CI=2.9000000953674316\n",
            "loading 5602...no preop cr\n",
            "loading 3555...1.01 -> 0.9, no AKI, min CI=1.600000023841858, max CI=3.9000000953674316\n",
            "loading 5607...0.88 -> 1.13, no AKI, min CI=0.8999999761581421, max CI=3.9000000953674316\n",
            "loading 1512...6.46 -> 5.07, no AKI, min CI=2.4000000953674316, max CI=3.9000000953674316\n",
            "loading 3562...0.75 -> 0.71, no AKI, min CI=0.800000011920929, max CI=3.5\n",
            "loading 3565...no postop cr\n",
            "loading 5613...0.73 -> 1.9, AKI, min CI=2.200000047683716, max CI=3.9000000953674316\n",
            "loading 3566...0.71 -> 0.66, no AKI, min CI=0.699999988079071, max CI=3.0999999046325684\n",
            "loading 1520...0.69 -> 0.64, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 5617...0.66 -> 0.62, no AKI, min CI=1.399999976158142, max CI=3.9000000953674316\n",
            "loading 5614...0.82 -> 0.87, no AKI, min CI=2.200000047683716, max CI=3.299999952316284\n",
            "loading 1521...0.85 -> 0.76, no AKI, min CI=3.4000000953674316, max CI=3.9000000953674316\n",
            "loading 1523...0.95 -> 0.91, no AKI, min CI=1.2999999523162842, max CI=3.0\n",
            "loading 3571...1.01 -> 1.14, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 3572...0.74 -> 0.67, no AKI, min CI=1.399999976158142, max CI=3.9000000953674316\n",
            "loading 3573...1.03 -> 0.87, no AKI, min CI=2.4000000953674316, max CI=3.9000000953674316\n",
            "loading 5624...0.78 -> 0.58, no AKI, min CI=2.200000047683716, max CI=3.9000000953674316\n",
            "loading 5621...3.93 -> 3.24, no AKI, min CI=0.699999988079071, max CI=3.9000000953674316\n",
            "loading 3576...0.89 -> 0.86, no AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 3578...0.95 -> 0.9, no AKI, min CI=2.0, max CI=3.4000000953674316\n",
            "loading 5626...9.26 -> 5.64, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 1536...no postop cr\n",
            "loading 5634...1.13 -> 1.24, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 5635...1.29 -> 1.21, no AKI, min CI=1.600000023841858, max CI=3.0999999046325684\n",
            "loading 1542...0.65 -> 0.59, no AKI, min CI=1.899999976158142, max CI=3.4000000953674316\n",
            "loading 1553...no preop cr\n",
            "loading 3602...no preop cr\n",
            "loading 1558...1.12 -> 1.05, no AKI, min CI=0.699999988079071, max CI=3.9000000953674316\n",
            "loading 5657...0.77 -> 0.55, no AKI, min CI=1.2999999523162842, max CI=2.200000047683716\n",
            "loading 1562...no preop cr\n",
            "loading 1564...0.74 -> 0.83, no AKI, min CI=2.799999952316284, max CI=3.9000000953674316\n",
            "loading 1566...0.67 -> 0.79, no AKI, min CI=1.2999999523162842, max CI=2.9000000953674316\n",
            "loading 3614...1.05 -> 1.03, no AKI, min CI=1.2000000476837158, max CI=2.5\n",
            "loading 3618...no preop cr\n",
            "loading 5669...18.1 -> 10.1, no AKI, min CI=2.5999999046325684, max CI=3.9000000953674316\n",
            "loading 1574...1.01 -> 0.96, no AKI, min CI=0.8999999761581421, max CI=3.299999952316284\n",
            "loading 5670...0.92 -> 0.99, no AKI, min CI=0.10000000149011612, max CI=3.9000000953674316\n",
            "loading 5671...0.86 -> 0.77, no AKI, min CI=1.399999976158142, max CI=3.9000000953674316\n",
            "loading 5675...no preop cr\n",
            "loading 3631...0.82 -> 0.66, no AKI, min CI=1.0, max CI=3.9000000953674316\n",
            "loading 5682...0.89 -> 1.0, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 5684...1.1 -> 0.91, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 1590...1.13 -> 2.78, AKI, min CI=3.0999999046325684, max CI=3.9000000953674316\n",
            "loading 5687...0.67 -> 0.66, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 5690...no postop cr\n",
            "loading 1597...0.85 -> 1.06, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 5694...no preop cr\n",
            "loading 5693...0.77 -> 0.95, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 5698...0.82 -> 1.46, AKI, min CI=2.4000000953674316, max CI=3.9000000953674316\n",
            "loading 3652...0.83 -> 0.9, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 3658...no preop cr\n",
            "loading 1614...0.61 -> 0.67, no AKI, min CI=1.399999976158142, max CI=2.9000000953674316\n",
            "loading 1615...0.96 -> 1.0, no AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 3663...no preop cr\n",
            "loading 5715...0.61 -> 0.62, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 5722...no preop cr\n",
            "loading 5724...0.57 -> 0.45, no AKI, min CI=1.5, max CI=3.0999999046325684\n",
            "loading 1630...0.73 -> 0.58, no AKI, min CI=1.5, max CI=2.799999952316284\n",
            "loading 1633...no postop cr\n",
            "loading 5729...0.65 -> 0.52, no AKI, min CI=1.399999976158142, max CI=2.700000047683716\n",
            "loading 3686...0.85 -> 0.74, no AKI, min CI=1.7000000476837158, max CI=3.0\n",
            "loading 1641...1.05 -> 1.89, AKI, min CI=1.399999976158142, max CI=2.700000047683716\n",
            "loading 1642...0.37 -> 0.44, no AKI, min CI=1.600000023841858, max CI=3.5\n",
            "loading 3694...0.65 -> 1.2, AKI, min CI=1.399999976158142, max CI=3.9000000953674316\n",
            "loading 1648...no preop cr\n",
            "loading 1652...1.02 -> 1.19, no AKI, min CI=0.8999999761581421, max CI=3.9000000953674316\n",
            "loading 3700...0.78 -> 1.52, AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 5750...5.17 -> 4.61, no AKI, min CI=1.399999976158142, max CI=3.9000000953674316\n",
            "loading 1656...no postop cr\n",
            "loading 1657...0.87 -> 0.88, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 3706...0.63 -> 0.65, no AKI, min CI=1.5, max CI=3.5\n",
            "loading 3710...7.18 -> 6.51, no AKI, min CI=1.100000023841858, max CI=3.9000000953674316\n",
            "loading 3719...0.76 -> 1.0, no AKI, min CI=1.2000000476837158, max CI=3.9000000953674316\n",
            "loading 1672...0.83 -> 1.58, AKI, min CI=1.399999976158142, max CI=3.9000000953674316\n",
            "loading 1673...1.24 -> 1.24, no AKI, min CI=2.5, max CI=3.9000000953674316\n",
            "loading 1674...0.86 -> 0.58, no AKI, min CI=1.2999999523162842, max CI=3.9000000953674316\n",
            "loading 1676...0.72 -> 0.6, no AKI, min CI=2.200000047683716, max CI=3.9000000953674316\n",
            "loading 3724...0.74 -> 0.62, no AKI, min CI=3.5, max CI=3.9000000953674316\n",
            "loading 5772...0.64 -> 0.68, no AKI, min CI=1.2000000476837158, max CI=3.0\n",
            "loading 5774...0.82 -> 0.81, no AKI, min CI=1.0, max CI=3.700000047683716\n",
            "loading 1681...0.86 -> 0.89, no AKI, min CI=1.7999999523162842, max CI=3.5999999046325684\n",
            "loading 5778...0.78 -> 0.86, no AKI, min CI=2.0999999046325684, max CI=3.0999999046325684\n",
            "loading 5777...0.53 -> 0.52, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 5780...0.84 -> 0.89, no AKI, min CI=2.200000047683716, max CI=3.9000000953674316\n",
            "loading 1688...1.09 -> 1.13, no AKI, min CI=1.2999999523162842, max CI=3.299999952316284\n",
            "loading 5787...no postop cr\n",
            "loading 3744...0.93 -> 1.19, no AKI, min CI=0.800000011920929, max CI=3.9000000953674316\n",
            "loading 1700...0.88 -> 0.86, no AKI, min CI=0.4000000059604645, max CI=3.9000000953674316\n",
            "loading 3752...0.73 -> 0.59, no AKI, min CI=2.200000047683716, max CI=3.9000000953674316\n",
            "loading 3753...0.99 -> 0.79, no AKI, min CI=0.800000011920929, max CI=3.9000000953674316\n",
            "loading 1706...0.93 -> 0.88, no AKI, min CI=2.0, max CI=3.4000000953674316\n",
            "loading 1708...1.17 -> 1.38, no AKI, min CI=0.800000011920929, max CI=2.5999999046325684\n",
            "loading 3757...0.61 -> 0.69, no AKI, min CI=2.299999952316284, max CI=3.200000047683716\n",
            "loading 1710...13.75 -> 12.66, no AKI, min CI=0.699999988079071, max CI=3.9000000953674316\n",
            "loading 5808...0.77 -> 0.64, no AKI, min CI=2.0999999046325684, max CI=2.9000000953674316\n",
            "loading 3761...no postop cr\n",
            "loading 5809...0.74 -> 0.55, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 3763...no preop cr\n",
            "loading 1716...8.88 -> 6.61, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 5816...12.89 -> 10.83, no AKI, min CI=2.200000047683716, max CI=3.9000000953674316\n",
            "loading 1726...0.68 -> 0.63, no AKI, min CI=2.799999952316284, max CI=3.9000000953674316\n",
            "loading 5825...1.03 -> 1.09, no AKI, min CI=0.800000011920929, max CI=2.799999952316284\n",
            "loading 1730...0.95 -> 1.71, AKI, min CI=0.4000000059604645, max CI=3.700000047683716\n",
            "loading 5831...0.6 -> 0.6, no AKI, min CI=1.7999999523162842, max CI=2.9000000953674316\n",
            "loading 1735...1.01 -> 1.09, no AKI, min CI=0.4000000059604645, max CI=3.799999952316284\n",
            "loading 3785...0.92 -> 0.78, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 5837...no ci\n",
            "loading 3791...0.94 -> 0.93, no AKI, min CI=0.8999999761581421, max CI=3.9000000953674316\n",
            "loading 5840...0.78 -> 0.55, no AKI, min CI=1.7000000476837158, max CI=2.5\n",
            "loading 5842...0.84 -> 1.06, no AKI, min CI=1.2000000476837158, max CI=3.9000000953674316\n",
            "loading 1749...0.85 -> 0.83, no AKI, min CI=2.5999999046325684, max CI=3.9000000953674316\n",
            "loading 3799...0.89 -> 1.12, no AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 5848...1.53 -> 1.27, no AKI, min CI=1.2999999523162842, max CI=3.299999952316284\n",
            "loading 1753...0.86 -> 0.49, no AKI, min CI=2.700000047683716, max CI=3.799999952316284\n",
            "loading 1754...0.58 -> 0.66, no AKI, min CI=2.0, max CI=2.9000000953674316\n",
            "loading 5853...no postop cr\n",
            "loading 5859...0.74 -> 0.68, no AKI, min CI=0.800000011920929, max CI=3.9000000953674316\n",
            "loading 1765...no postop cr\n",
            "loading 3814...0.55 -> 0.52, no AKI, min CI=2.0, max CI=3.799999952316284\n",
            "loading 1771...1.15 -> 1.0, no AKI, min CI=0.800000011920929, max CI=3.5\n",
            "loading 3819...1.1 -> 0.85, no AKI, min CI=0.800000011920929, max CI=2.700000047683716\n",
            "loading 1773...0.64 -> 0.62, no AKI, min CI=2.4000000953674316, max CI=3.799999952316284\n",
            "loading 5870...1.22 -> 1.3, no AKI, min CI=1.2999999523162842, max CI=2.9000000953674316\n",
            "loading 5872...0.74 -> 1.0, no AKI, min CI=2.5, max CI=3.9000000953674316\n",
            "loading 3825...0.98 -> 0.85, no AKI, min CI=0.30000001192092896, max CI=3.5999999046325684\n",
            "loading 3828...no preop cr\n",
            "loading 1784...0.84 -> 0.88, no AKI, min CI=2.5999999046325684, max CI=3.9000000953674316\n",
            "loading 1785...0.65 -> 1.21, AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 5882...0.98 -> 0.9, no AKI, min CI=3.299999952316284, max CI=3.9000000953674316\n",
            "loading 5884...3.24 -> 2.98, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 5888...0.59 -> 0.65, no AKI, min CI=1.2000000476837158, max CI=2.299999952316284\n",
            "loading 3842...0.7 -> 0.68, no AKI, min CI=1.600000023841858, max CI=3.200000047683716\n",
            "loading 3844...0.89 -> 0.85, no AKI, min CI=1.7000000476837158, max CI=3.4000000953674316\n",
            "loading 3846...1.05 -> 0.9, no AKI, min CI=1.7999999523162842, max CI=3.799999952316284\n",
            "loading 3849...1.15 -> 1.06, no AKI, min CI=1.2000000476837158, max CI=3.9000000953674316\n",
            "loading 3850...0.98 -> 0.97, no AKI, min CI=1.2999999523162842, max CI=3.9000000953674316\n",
            "loading 1803...0.62 -> 0.6, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 1802...0.77 -> 0.62, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 5898...8.68 -> 7.07, no AKI, min CI=2.5999999046325684, max CI=3.9000000953674316\n",
            "loading 1804...2.21 -> 1.79, no AKI, min CI=2.4000000953674316, max CI=3.9000000953674316\n",
            "loading 1805...0.79 -> 0.54, no AKI, min CI=1.5, max CI=2.799999952316284\n",
            "loading 3863...0.91 -> 0.87, no AKI, min CI=0.6000000238418579, max CI=3.4000000953674316\n",
            "loading 5912...0.8 -> 0.82, no AKI, min CI=0.6000000238418579, max CI=2.799999952316284\n",
            "loading 5914...no preop cr\n",
            "loading 1820...0.91 -> 1.04, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 5917...0.34 -> 0.44, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 3870...0.6 -> 0.65, no AKI, min CI=2.4000000953674316, max CI=3.9000000953674316\n",
            "loading 5921...no postop cr\n",
            "loading 1826...0.97 -> 0.86, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 3878...1.16 -> 1.13, no AKI, min CI=1.7999999523162842, max CI=3.799999952316284\n",
            "loading 3879...0.87 -> 0.9, no AKI, min CI=1.7999999523162842, max CI=3.5999999046325684\n",
            "loading 5928...0.74 -> 0.73, no AKI, min CI=1.899999976158142, max CI=3.299999952316284\n",
            "loading 1836...0.93 -> 0.9, no AKI, min CI=1.100000023841858, max CI=3.9000000953674316\n",
            "loading 5936...no preop cr\n",
            "loading 5937...no preop cr\n",
            "loading 3889...1.07 -> 0.89, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 3890...7.53 -> 8.04, no AKI, min CI=1.2000000476837158, max CI=3.9000000953674316\n",
            "loading 1844...0.8 -> 0.76, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 5940...1.28 -> 1.23, no AKI, min CI=1.899999976158142, max CI=3.799999952316284\n",
            "loading 1846...0.95 -> 1.01, no AKI, min CI=1.399999976158142, max CI=2.799999952316284\n",
            "loading 3895...no preop cr\n",
            "loading 5946...1.1 -> 1.12, no AKI, min CI=1.399999976158142, max CI=2.700000047683716\n",
            "loading 5948...0.89 -> 0.83, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 1852...0.73 -> 0.84, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 3902...0.84 -> 0.79, no AKI, min CI=1.600000023841858, max CI=3.9000000953674316\n",
            "loading 5951...0.87 -> 0.72, no AKI, min CI=2.299999952316284, max CI=3.700000047683716\n",
            "loading 1855...1.52 -> 1.51, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 3904...0.72 -> 0.7, no AKI, min CI=1.2999999523162842, max CI=3.9000000953674316\n",
            "loading 1865...no postop cr\n",
            "loading 1866...no postop cr\n",
            "loading 5964...0.81 -> 0.79, no AKI, min CI=0.800000011920929, max CI=3.9000000953674316\n",
            "loading 5966...0.85 -> 0.84, no AKI, min CI=1.2999999523162842, max CI=3.0999999046325684\n",
            "loading 5971...0.53 -> 0.55, no AKI, min CI=1.2999999523162842, max CI=3.799999952316284\n",
            "loading 3925...1.56 -> 1.32, no AKI, min CI=0.8999999761581421, max CI=3.0\n",
            "loading 5973...0.71 -> 0.55, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 5977...no preop cr\n",
            "loading 3930...0.6 -> 1.44, AKI, min CI=2.5, max CI=3.9000000953674316\n",
            "loading 3935...1.06 -> 0.99, no AKI, min CI=1.7000000476837158, max CI=3.4000000953674316\n",
            "loading 5983...0.54 -> 1.2, AKI, min CI=1.2999999523162842, max CI=3.9000000953674316\n",
            "loading 1891...0.85 -> 0.71, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 1896...5.46 -> 5.42, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 1897...1.08 -> 1.1, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 1898...0.69 -> 0.66, no AKI, min CI=2.200000047683716, max CI=3.0999999046325684\n",
            "loading 1900...0.6 -> 0.45, no AKI, min CI=0.20000000298023224, max CI=3.9000000953674316\n",
            "loading 3950...0.66 -> 0.58, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 1903...0.79 -> 0.71, no AKI, min CI=1.0, max CI=3.9000000953674316\n",
            "loading 1912...0.81 -> 0.67, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 6009...0.9 -> 0.96, no AKI, min CI=0.8999999761581421, max CI=3.9000000953674316\n",
            "loading 3962...0.98 -> 0.92, no AKI, min CI=2.0999999046325684, max CI=3.799999952316284\n",
            "loading 6015...0.81 -> 0.76, no AKI, min CI=0.8999999761581421, max CI=3.9000000953674316\n",
            "loading 3967...0.82 -> 0.62, no AKI, min CI=1.2000000476837158, max CI=2.4000000953674316\n",
            "loading 6017...0.6 -> 0.71, no AKI, min CI=2.799999952316284, max CI=3.9000000953674316\n",
            "loading 3971...0.59 -> 0.53, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 3972...0.75 -> 0.68, no AKI, min CI=0.699999988079071, max CI=2.5999999046325684\n",
            "loading 3973...0.5 -> 0.49, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 6022...0.68 -> 0.73, no AKI, min CI=0.8999999761581421, max CI=2.700000047683716\n",
            "loading 3975...1.04 -> 0.92, no AKI, min CI=1.2999999523162842, max CI=3.0\n",
            "loading 6027...0.97 -> 0.68, no AKI, min CI=0.800000011920929, max CI=3.200000047683716\n",
            "loading 6029...no postop cr\n",
            "loading 1935...0.76 -> 0.74, no AKI, min CI=1.600000023841858, max CI=2.700000047683716\n",
            "loading 1936...0.7 -> 0.96, no AKI, min CI=1.0, max CI=3.700000047683716\n",
            "loading 6031...0.78 -> 0.72, no AKI, min CI=1.399999976158142, max CI=3.9000000953674316\n",
            "loading 6032...1.31 -> 1.3, no AKI, min CI=2.5999999046325684, max CI=3.9000000953674316\n",
            "loading 1938...0.64 -> 0.65, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 1941...0.61 -> 0.75, no AKI, min CI=1.899999976158142, max CI=3.9000000953674316\n",
            "loading 1942...0.73 -> 0.71, no AKI, min CI=1.2999999523162842, max CI=2.5\n",
            "loading 6037...4.91 -> 3.63, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 3991...0.77 -> 0.84, no AKI, min CI=2.299999952316284, max CI=3.9000000953674316\n",
            "loading 3994...0.83 -> 0.88, no AKI, min CI=0.8999999761581421, max CI=3.9000000953674316\n",
            "loading 6042...0.84 -> 0.79, no AKI, min CI=1.2999999523162842, max CI=3.9000000953674316\n",
            "loading 1949...no preop cr\n",
            "loading 1957...2.62 -> 2.42, no AKI, min CI=2.700000047683716, max CI=3.9000000953674316\n",
            "loading 1959...1.0 -> 0.99, no AKI, min CI=0.6000000238418579, max CI=3.5\n",
            "loading 6059...0.8 -> 0.85, no AKI, min CI=1.5, max CI=3.299999952316284\n",
            "loading 4013...no preop cr\n",
            "loading 1969...0.74 -> 0.76, no AKI, min CI=0.800000011920929, max CI=3.0\n",
            "loading 4020...0.67 -> 0.67, no AKI, min CI=1.600000023841858, max CI=3.799999952316284\n",
            "loading 1973...1.28 -> 1.48, no AKI, min CI=2.4000000953674316, max CI=3.9000000953674316\n",
            "loading 4021...0.89 -> 0.87, no AKI, min CI=1.899999976158142, max CI=3.5\n",
            "loading 6071...0.81 -> 0.68, no AKI, min CI=1.600000023841858, max CI=3.9000000953674316\n",
            "loading 4024...no postop cr\n",
            "loading 4022...0.46 -> 0.43, no AKI, min CI=2.700000047683716, max CI=3.9000000953674316\n",
            "loading 6074...1.1 -> 1.26, no AKI, min CI=1.2000000476837158, max CI=3.9000000953674316\n",
            "loading 4027...0.91 -> 0.9, no AKI, min CI=2.5999999046325684, max CI=3.9000000953674316\n",
            "loading 6077...no preop cr\n",
            "loading 4033...7.83 -> 7.3, no AKI, min CI=1.7000000476837158, max CI=3.9000000953674316\n",
            "loading 4034...1.57 -> 2.67, AKI, min CI=3.200000047683716, max CI=3.9000000953674316\n",
            "loading 1988...0.7 -> 0.68, no AKI, min CI=2.4000000953674316, max CI=3.9000000953674316\n",
            "loading 4036...0.84 -> 0.9, no AKI, min CI=1.5, max CI=3.299999952316284\n",
            "loading 4040...1.05 -> 0.94, no AKI, min CI=2.0999999046325684, max CI=3.5999999046325684\n",
            "loading 6089...no postop cr\n",
            "loading 4042...1.01 -> 1.04, no AKI, min CI=1.2999999523162842, max CI=3.5999999046325684\n",
            "loading 4048...0.85 -> 0.69, no AKI, min CI=0.8999999761581421, max CI=3.299999952316284\n",
            "loading 2000...1.29 -> 1.86, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 6104...0.77 -> 0.77, no AKI, min CI=2.0999999046325684, max CI=3.9000000953674316\n",
            "loading 2010...0.83 -> 0.99, no AKI, min CI=2.0, max CI=3.9000000953674316\n",
            "loading 4060...0.68 -> 0.61, no AKI, min CI=0.8999999761581421, max CI=3.0\n",
            "loading 2014...0.65 -> 0.69, no AKI, min CI=1.2000000476837158, max CI=3.9000000953674316\n",
            "loading 2016...0.85 -> 0.92, no AKI, min CI=1.7999999523162842, max CI=3.9000000953674316\n",
            "loading 2017...0.56 -> 0.51, no AKI, min CI=1.100000023841858, max CI=2.9000000953674316\n",
            "loading 6114...0.78 -> 0.81, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 6119...no preop cr\n",
            "loading 6121...0.91 -> 0.72, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "loading 2028...0.66 -> 0.55, no AKI, min CI=1.100000023841858, max CI=3.799999952316284\n",
            "loading 6126...0.68 -> 0.67, no AKI, min CI=2.200000047683716, max CI=3.9000000953674316\n",
            "loading 6133...1.52 -> 1.62, no AKI, min CI=1.100000023841858, max CI=3.9000000953674316\n",
            "loading 6135...0.64 -> 0.56, no AKI, min CI=2.5999999046325684, max CI=3.299999952316284\n",
            "loading 2040...no preop cr\n",
            "loading 2044...0.58 -> 0.56, no AKI, min CI=1.5, max CI=3.9000000953674316\n",
            "35 AKI/ 808 cases 4.3%\n"
          ]
        }
      ],
      "source": [
        "# Set blood pressure threshold\n",
        "thresholds = np.arange(0, 4, 0.1)\n",
        "\n",
        "# Save the final result\n",
        "rows = []\n",
        "for caseid in caseids:\n",
        "    print('loading {}...'.format(caseid), flush=True, end='')\n",
        "\n",
        "    # Column ['anend'] : anesthesia end time\n",
        "    aneend = df_cases[(df_cases['caseid'] == caseid)]['aneend'].values[0]\n",
        "\n",
        "    # Last creatinine concentration before surgery\n",
        "    preop_cr = df_labs[(df_labs['caseid'] == caseid) & (df_labs['dt'] < 0) & (df_labs['name'] == 'cr')].sort_values(by=['dt'], axis=0, ascending=False)['result'].values.flatten()\n",
        "    if len(preop_cr) == 0:\n",
        "        print('no preop cr')\n",
        "        continue\n",
        "    preop_cr = preop_cr[0]\n",
        "\n",
        "    # Maximum creatinine concentration within 48 hours after surgery\n",
        "    postop_cr = df_labs[(df_labs['caseid'] == caseid) & (df_labs['dt'] > aneend) &\n",
        "        (df_labs['dt'] < aneend + 48 * 3600) & (df_labs['name'] == 'cr')]['result'].max(skipna=True)\n",
        "    if not postop_cr or np.isnan(postop_cr):\n",
        "        print('no postop cr')\n",
        "        continue\n",
        "\n",
        "    # KDIGO stage I\n",
        "    aki = postop_cr > preop_cr * 1.5\n",
        "\n",
        "    # Blood pressure during surgery\n",
        "    cis = vitaldb.load_case(caseid, '/CI').flatten()\n",
        "    cis = cis[~np.isnan(cis)]\n",
        "    cis = cis[(cis > 0.1) & (cis < 4)]\n",
        "    if len(cis) < 10:\n",
        "        print('no ci')\n",
        "        continue\n",
        "\n",
        "    # Calculate the percentage that stays for the time as increasing the blood pressure by 1 unit.\n",
        "    row = {'aki':aki}\n",
        "    for threshold in thresholds:\n",
        "        row[f'under{threshold}'] = np.nanmean(cis < threshold) * 100\n",
        "\n",
        "    # Append the result into row\n",
        "    rows.append(row)\n",
        "\n",
        "    print(f'{preop_cr} -> {postop_cr}, {\"AKI\" if aki else \"no AKI\"}, min CI={min(cis)}, max CI={max(cis)}')\n",
        "\n",
        "df = pd.DataFrame(rows)\n",
        "print(f'{df[\"aki\"].sum()} AKI/ {len(df)} cases {df[\"aki\"].mean() * 100:.1f}%')"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "YYHU8na4WV-0"
      },
      "source": [
        "## Calculate the odds ratio for AKI prediction depending on blood pressure threshold\n",
        "\n",
        "- Let's calculate the odds ratio regarding how much each BP class increases the risk of AKI, using the above-measured values."
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "print(thresholds)"
      ],
      "metadata": {
        "id": "vxGUAaSY7R5p",
        "outputId": "64363b8d-7a1c-4a7f-dd7b-5b4b08c87b25",
        "colab": {
          "base_uri": "https://localhost:8080/"
        }
      },
      "execution_count": 21,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "[0.  0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.  1.1 1.2 1.3 1.4 1.5 1.6 1.7\n",
            " 1.8 1.9 2.  2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3.  3.1 3.2 3.3 3.4 3.5\n",
            " 3.6 3.7 3.8 3.9]\n"
          ]
        }
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 26,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "X5O0yxHJWcjh",
        "outputId": "38823016-e7e3-4a40-a23b-ef14d984f740"
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "Optimization terminated successfully.\n",
            "         Current function value: 0.178346\n",
            "         Iterations 7\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.178346\n",
            "         Iterations 7\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.175587\n",
            "         Iterations 7\n",
            "under0.2\tb=0.248, exp(b)=1.281, pval=0.018\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.174874\n",
            "         Iterations 7\n",
            "under0.30000000000000004\tb=0.149, exp(b)=1.161, pval=0.036\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.174802\n",
            "         Iterations 7\n",
            "under0.4\tb=0.135, exp(b)=1.144, pval=0.030\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.174555\n",
            "         Iterations 7\n",
            "under0.5\tb=0.129, exp(b)=1.137, pval=0.026\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.174458\n",
            "         Iterations 7\n",
            "under0.6000000000000001\tb=0.123, exp(b)=1.131, pval=0.022\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.174464\n",
            "         Iterations 7\n",
            "under0.7000000000000001\tb=0.109, exp(b)=1.116, pval=0.016\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.174464\n",
            "         Iterations 7\n",
            "under0.8\tb=0.109, exp(b)=1.116, pval=0.016\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.174146\n",
            "         Iterations 7\n",
            "under0.9\tb=0.101, exp(b)=1.106, pval=0.008\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.174146\n",
            "         Iterations 7\n",
            "under1.0\tb=0.101, exp(b)=1.106, pval=0.008\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.174171\n",
            "         Iterations 7\n",
            "under1.1\tb=0.094, exp(b)=1.099, pval=0.006\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.174250\n",
            "         Iterations 7\n",
            "under1.2000000000000002\tb=0.084, exp(b)=1.088, pval=0.005\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.175993\n",
            "         Iterations 7\n",
            "under1.3\tb=0.049, exp(b)=1.050, pval=0.023\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.177066\n",
            "         Iterations 7\n",
            "under1.4000000000000001\tb=0.031, exp(b)=1.031, pval=0.100\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.177066\n",
            "         Iterations 7\n",
            "under1.5\tb=0.031, exp(b)=1.031, pval=0.100\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.177748\n",
            "         Iterations 7\n",
            "under1.6\tb=0.018, exp(b)=1.018, pval=0.278\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.178175\n",
            "         Iterations 7\n",
            "under1.7000000000000002\tb=0.008, exp(b)=1.008, pval=0.580\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.178304\n",
            "         Iterations 7\n",
            "under1.8\tb=-0.003, exp(b)=0.997, pval=0.801\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.178055\n",
            "         Iterations 7\n",
            "under1.9000000000000001\tb=-0.007, exp(b)=0.993, pval=0.514\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.178055\n",
            "         Iterations 7\n",
            "under2.0\tb=-0.007, exp(b)=0.993, pval=0.514\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.177143\n",
            "         Iterations 7\n",
            "under2.1\tb=-0.010, exp(b)=0.990, pval=0.197\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.177143\n",
            "         Iterations 7\n",
            "under2.2\tb=-0.010, exp(b)=0.990, pval=0.197\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.176715\n",
            "         Iterations 7\n",
            "under2.3000000000000003\tb=-0.010, exp(b)=0.990, pval=0.127\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.176715\n",
            "         Iterations 7\n",
            "under2.4000000000000004\tb=-0.010, exp(b)=0.990, pval=0.127\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.177044\n",
            "         Iterations 7\n",
            "under2.5\tb=-0.008, exp(b)=0.992, pval=0.164\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.177011\n",
            "         Iterations 7\n",
            "under2.6\tb=-0.007, exp(b)=0.993, pval=0.151\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.177011\n",
            "         Iterations 7\n",
            "under2.7\tb=-0.007, exp(b)=0.993, pval=0.151\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.176956\n",
            "         Iterations 7\n",
            "under2.8000000000000003\tb=-0.007, exp(b)=0.993, pval=0.137\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.176956\n",
            "         Iterations 7\n",
            "under2.9000000000000004\tb=-0.007, exp(b)=0.993, pval=0.137\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.177171\n",
            "         Iterations 7\n",
            "under3.0\tb=-0.007, exp(b)=0.993, pval=0.168\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.177226\n",
            "         Iterations 7\n",
            "under3.1\tb=-0.007, exp(b)=0.993, pval=0.173\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.177226\n",
            "         Iterations 7\n",
            "under3.2\tb=-0.007, exp(b)=0.993, pval=0.173\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.177161\n",
            "         Iterations 7\n",
            "under3.3000000000000003\tb=-0.008, exp(b)=0.992, pval=0.155\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.177161\n",
            "         Iterations 7\n",
            "under3.4000000000000004\tb=-0.008, exp(b)=0.992, pval=0.155\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.177102\n",
            "         Iterations 7\n",
            "under3.5\tb=-0.009, exp(b)=0.991, pval=0.141\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.176958\n",
            "         Iterations 7\n",
            "under3.6\tb=-0.012, exp(b)=0.988, pval=0.112\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.176958\n",
            "         Iterations 7\n",
            "under3.7\tb=-0.012, exp(b)=0.988, pval=0.112\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.177698\n",
            "         Iterations 7\n",
            "under3.8000000000000003\tb=-0.020, exp(b)=0.980, pval=0.268\n",
            "Optimization terminated successfully.\n",
            "         Current function value: 0.177698\n",
            "         Iterations 7\n",
            "under3.9000000000000004\tb=-0.020, exp(b)=0.980, pval=0.268\n"
          ]
        }
      ],
      "source": [
        "import statsmodels.api as sm\n",
        "\n",
        "# Get odd ration using univariate logistic regression\n",
        "df['intercept'] = 1\n",
        "df['aki'] = df['aki'].astype(bool)\n",
        "odd_ratios = []\n",
        "for threshold in thresholds:\n",
        "    c = f'under{threshold}'\n",
        "    try:\n",
        "        model = sm.Logit(df['aki'], df[['intercept', c]])\n",
        "        res = model.fit()\n",
        "    except:\n",
        "        odd_ratios.append(np.nan)\n",
        "        continue\n",
        "    b = res.params[c]\n",
        "    pval = res.pvalues[c]\n",
        "    odd_ratios.append(np.exp(b))\n",
        "    print(f'{c}\\tb={b:.3f}, exp(b)={np.exp(b):.3f}, pval={pval:.3f}')"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "baARXo6kWeO4"
      },
      "source": [
        "## Draw a graph\n",
        "- We can find out that the odds ratio of AKI increases as the time that MBP stays below 60-65mmHg increases.\n"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 27,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 430
        },
        "id": "mut5q3MeWnKQ",
        "outputId": "7fd5d260-7cc2-45f6-edcc-14bbe4c1bbd7"
      },
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "text/plain": [
              "<Figure size 640x480 with 1 Axes>"
            ],
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAi0AAAGdCAYAAADey0OaAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA2zUlEQVR4nO3deXhU5d3/8c+ZLJMASSAkZJGwiYBsYZMYfRRQMCJPKrYVRKupLa22aEHaUvjVitY+YlsVFyhqrabaKoIKVkEQUaAIqARGQBZZAgRIwp4Vss35/RESjRCYCTNzZibv13WdC3Nyn+F7OMT5cOa+v8cwTdMUAACAn7NZXQAAAIArCC0AACAgEFoAAEBAILQAAICAQGgBAAABgdACAAACAqEFAAAEBEILAAAICKFWF+ApTqdThw4dUlRUlAzDsLocAADgAtM0VVJSouTkZNls57+XEjSh5dChQ0pJSbG6DAAA0AR5eXlq3779eccETWiJioqSVHvS0dHRFlcDAABcUVxcrJSUlPr38fMJmtBS95FQdHQ0oQUAgADjytQOJuICAICAQGgBAAABgdACAAACAqEFAAAEBEILAAAICIQWAAAQEAgtAAAgIBBaAABAQCC0AACAgEBoAQAAAYHQAgAAAgKhBQAABARCywU8/dHXmvbOJh0vq7S6FAAAmjVCywX8+7P9euPzPB06ecrqUgAAaNYILRcQ18ouSTpaWmFxJQAANG+ElguIaxUuSTpaysdDAABYidByAfFn7rQcKeFOCwAAViK0XEB8FB8PAQDgDwgtF8CcFgAA/AOh5QLiourmtBBaAACwEqHlAuJbRUhiTgsAAFYjtFzAN3daWD0EAICVCC0XUDen5UR5papqnBZXAwBA80VouYA2LcIVYjNkmqKVPwAAFiK0XECIzVBsy9qPiJjXAgCAdQgtLmDZMwAA1iO0uIBW/gAAWI/Q4oK6rrh8PAQAgHUILS6I5+MhAAAsR2hxAXNaAACwHqHFBXw8BACA9QgtLuBOCwAA1iO0uIBW/gAAWI/Q4gJa+QMAYD1Ciwto5Q8AgPUILS6glT8AANYjtLiIybgAAFiL0OIilj0DAGAtQouLeP4QAADWIrS4iFb+AABYi9Dioro5LXw8BACANQgtLqqb08KdFgAArEFocRGrhwAAsBahxUW08gcAwFpuh5ZVq1YpMzNTycnJMgxDCxcuPO/4d955RyNGjFB8fLyio6OVnp6upUuXNhjz8MMPyzCMBluPHj3cLc2r6ibiHi+jlT8AAFZwO7SUlZUpNTVVs2fPdmn8qlWrNGLECC1evFg5OTkaNmyYMjMztXHjxgbjevXqpfz8/Ppt9erV7pbmVXWt/CVa+QMAYIVQdw8YOXKkRo4c6fL4p59+usHXjz32mN59912999576t+//zeFhIYqMTHR3XJ8xnamlf+RkgodKalQQnSE1SUBANCs+HxOi9PpVElJiWJjYxvs37lzp5KTk9WlSxfdcccd2r9//3lfp6KiQsXFxQ02b6tf9sxkXAAAfM7noeWJJ55QaWmpxowZU78vLS1N2dnZWrJkiebMmaPc3Fxdc801KikpafR1ZsyYoZiYmPotJSXF67XXL3umVwsAAD7n09Dy+uuv65FHHtG8efPUrl27+v0jR47Urbfeqr59+yojI0OLFy/WyZMnNW/evEZfa9q0aSoqKqrf8vLyvF4/rfwBALCO23Nammru3LkaP3685s+fr+HDh593bOvWrdWtWzft2rWr0TF2u112u93TZZ4XrfwBALCOT+60vPHGG7r77rv1xhtvaNSoURccX1paqt27dyspKckH1bmOJz0DAGAdt++0lJaWNrgDkpubK4fDodjYWHXo0EHTpk3TwYMH9eqrr0qq/UgoKytLzzzzjNLS0lRQUCBJioyMVExMjCTpN7/5jTIzM9WxY0cdOnRI06dPV0hIiMaNG+eJc/QYuuICAGAdt++0rF+/Xv37969frjx58mT1799fDz30kCQpPz+/wcqfF198UdXV1ZowYYKSkpLqt4kTJ9aPOXDggMaNG6fu3btrzJgxatu2rdatW6f4+PiLPT+PIrQAAGAdwzRN0+oiPKG4uFgxMTEqKipSdHS0V36P7QXFuvHp/6pNizBtfOgGr/weAAA0J+68f/PsITfUTcQ9UV5FK38AAHyM0OIGWvkDAGAdQosb6lr5S6wgAgDA1wgtboqnlT8AAJYgtLgpjlb+AABYgtDiJlr5AwBgDUKLm+iKCwCANQgtbuL5QwAAWIPQ4ia64gIAYA1Ci5sILQAAWIPQ4ibmtAAAYA1Ci5vqVg/Ryh8AAN8itLiJVv4AAFiD0OImm81QW1r5AwDgc4SWJoijlT8AAD5HaGkCWvkDAOB7hJYmoJU/AAC+R2hpApY9AwDge4SWJqCVPwAAvkdoaQK64gIA4HuElibg4yEAAHyP0NIE3GkBAMD3CC1NQCt/AAB8j9DSBN9u5X+MZc8AAPgEoaUJvt3Kn4+IAADwDUJLE9HKHwAA3yK0NBGt/AEA8C1CSxPFc6cFAACfIrQ0UVzUmTktJUzEBQDAFwgtTUQrfwAAfIvQ0kT1E3GZ0wIAgE8QWpqorpU/d1oAAPANQksT0cofAADfIrQ0Ea38AQDwLUJLE9HKHwAA3yK0NBGt/AEA8C1Cy0WglT8AAL5DaLkIda38WfYMAID3EVouAg3mAADwHULLRaCVPwAAvkNouQjcaQEAwHcILRchnjktAAD4DKHlItAVFwAA3yG0XARCCwAAvkNouQh1Hw/Ryh8AAO9zO7SsWrVKmZmZSk5OlmEYWrhw4XnHv/POOxoxYoTi4+MVHR2t9PR0LV269Kxxs2fPVqdOnRQREaG0tDR9/vnn7pbmc60jw2jlDwCAj7gdWsrKypSamqrZs2e7NH7VqlUaMWKEFi9erJycHA0bNkyZmZnauHFj/Zg333xTkydP1vTp07VhwwalpqYqIyNDhw8fdrc8n6KVPwAAvmOYpmk2+WDD0IIFCzR69Gi3juvVq5fGjh2rhx56SJKUlpamK664QrNmzZIkOZ1OpaSk6P7779fUqVNdes3i4mLFxMSoqKhI0dHRbtVzMW565r/aml+sV+6+QsO6t/PZ7wsAQDBw5/3b53NanE6nSkpKFBsbK0mqrKxUTk6Ohg8f/k1RNpuGDx+utWvXNvo6FRUVKi4ubrBZgWXPAAD4hs9DyxNPPKHS0lKNGTNGknT06FHV1NQoISGhwbiEhAQVFBQ0+jozZsxQTExM/ZaSkuLVuhvDCiIAAHzDp6Hl9ddf1yOPPKJ58+apXbuL+yhl2rRpKioqqt/y8vI8VKV7aOUPAIBvhPrqN5o7d67Gjx+v+fPnN/goKC4uTiEhISosLGwwvrCwUImJiY2+nt1ul91u91q9rqpr5X+EOy0AAHiVT+60vPHGG7r77rv1xhtvaNSoUQ2+Fx4eroEDB2r58uX1+5xOp5YvX6709HRflHdR6ua0HGVOCwAAXuX2nZbS0lLt2rWr/uvc3Fw5HA7FxsaqQ4cOmjZtmg4ePKhXX31VUu1HQllZWXrmmWeUlpZWP08lMjJSMTExkqTJkycrKytLgwYN0uDBg/X000+rrKxMd999tyfO0auY0wIAgG+4HVrWr1+vYcOG1X89efJkSVJWVpays7OVn5+v/fv313//xRdfVHV1tSZMmKAJEybU768bL0ljx47VkSNH9NBDD6mgoED9+vXTkiVLzpqc648ILQAA+MZF9WnxJ1b1aTleVqkBjy6TJO38v5EKC+HJCAAAuMqv+7QEG1r5AwDgG4SWi0QrfwAAfIPQ4gF0xQUAwPsILR4QR68WAAC8jtDiAawgAgDA+wgtHlDXyp+PhwAA8B5CiwfE199pYfUQAADeQmjxAFr5AwDgfYQWD2BOCwAA3kdo8YD6Jc+EFgAAvIbQ4gF1d1pOllepqsZpcTUAAAQnQosH0MofAADvI7R4wLdb+bPsGQAA7yC0eEj9CiLmtQAA4BWEFg+hlT8AAN5FaPEQlj0DAOBdhBYP4UnPAAB4F6HFQ+Ja1U7EpZU/AADeQWjxEFr5AwDgXYQWD2EiLgAA3kVo8RCWPAMA4F2EFg+hlT8AAN5FaPEQWvkDAOBdhBYPsdmM+hVELHsGAMDzCC0eRIM5AAC8h9DiQawgAgDAewgtHkRXXAAAvIfQ4kF8PAQAgPcQWjyIVv4AAHgPocWDaOUPAID3EFo8KJ6JuAAAeA2hxYPiaOUPAIDXEFo86Nut/CuraeUPAIAnEVo8qHVkmELrWvmXcbcFAABPIrR4kM1mqG3dCqISVhABAOBJhBYPo1cLAADeQWjxsOTWkZKkD7bkW1wJAADBhdDiYT+/toskad76A/p011GLqwEAIHgQWjzsik6xuvPKjpKkqe9sUnlltcUVAQAQHAgtXjDlxu5KjolQ3vFTeurDr60uBwCAoEBo8YKoiDD93y19JEkvf5qrjftPWFwRAACBj9DiJcN6tNMt/S+R05Smvr2ZZnMAAFwkQosX/eF/eyq2Zbh2FJZozordVpcDAEBAI7R4UWzLcD38vV6SpFmf7NTXhSUWVwQAQOAitHhZZt8kDb+8napqTE15a5NqnKbVJQEAEJDcDi2rVq1SZmamkpOTZRiGFi5ceN7x+fn5uv3229WtWzfZbDZNmjTprDHZ2dkyDKPBFhER4W5pfskwDD06urei7KFy5J3UP9fstbokAAACktuhpaysTKmpqZo9e7ZL4ysqKhQfH68HH3xQqampjY6Ljo5Wfn5+/bZv3z53S/NbSTGRmnpTD0nSX5fuUN7xcosrAgAg8IS6e8DIkSM1cuRIl8d36tRJzzzzjCTp5ZdfbnScYRhKTEx0t5yAMe6KDvqP45A+yz2uae9s1ms/HSzDMKwuCwCAgOE3c1pKS0vVsWNHpaSk6Oabb9ZXX3113vEVFRUqLi5usPkzm83Q4z/oK3uoTat3HdVbOQesLgkAgIDiF6Gle/fuevnll/Xuu+/qX//6l5xOp6666iodOND4G/uMGTMUExNTv6WkpPiw4qbpHNdSk0d0kyQ9+v5WHS45bXFFAAAEDr8ILenp6brrrrvUr18/DRkyRO+8847i4+P1wgsvNHrMtGnTVFRUVL/l5eX5sOKm++n/dFafS2JUfLpa0989/90kAADwDb8ILd8VFham/v37a9euXY2Osdvtio6ObrAFgtAQm/78g74KtRn6YEuBlmzJt7okAAACgl+GlpqaGm3evFlJSUlWl+IVPZOjde+QSyVJf3j3KxWVV1lcEQAA/s/t0FJaWiqHwyGHwyFJys3NlcPh0P79+yXVfmxz1113NTimbnxpaamOHDkih8OhrVu31n//j3/8oz788EPt2bNHGzZs0I9+9CPt27dP48ePv4hT82/3XddVXeJb6khJhf5v8dYLHwAAQDPn9pLn9evXa9iwYfVfT548WZKUlZWl7Oxs5efn1weYOv3796//75ycHL3++uvq2LGj9u7dK0k6ceKEfvazn6mgoEBt2rTRwIEDtWbNGvXs2bMp5xQQIsJC9Jcf9NWtL6zVvPUHdEWnWP1wYHuWQQMA0AjDNM2g6CtfXFysmJgYFRUVBcz8Fkl6+D9fKftMl9xh3eP1p1v66JLWkdYWBQCAj7jz/u2Xc1qak9+PulyTR3RTeIhNn+w4ohFPrdTLq3N5RhEAAN9BaLFYWIhNv7r+Mi2eeI0Gd4pVeWWN/vj+Vn1/zhpty/fvhnkAAPgSocVPdG3XSnN/fqUeu6WPoiJC9WXeSWU+t1p/WbJdp6tqrC4PAADLEVr8iM1m6Pa0Dlo+eYhG9k5UtdPU31bs1o1Pr9Ka3UetLg8AAEsRWvxQu+gIzfnRQL1450AlRNu191i5bv/7Z5ry1pc6WV5pdXkAAFiC0OLHbuiVqGWTh+hHV3aQJM1bf0DDn1qp9748pCBZ9AUAgMsILX4uOiJMfxrdR2/dm66u7VrpaGml7n9jo36S/YXyjpdbXR4AAD5DaAkQgzrFatGv/kcPDP/W8uiZKzX7k12qrHZaXR4AAF5HaAkg9tAQTRxeuzw6vUtbna5y6q9Ld2jUs//V57nHrS4PAACvIrQEoK7tWun1n6Vp5thUtW0Zrp2HSzXmhbX67fwvdbyMiboAgOBEaAlQhmHolv7t9fGvh+r2tNqJuvNzDui6J1do3hd5ctJRFwAQZAgtAS6mRZgeu6WP3v7FVeqRGKWT5VWa8vYm3fbiOn1dWGJ1eQAAeAyhJUgM7NhG793/P/r9TZcrMixEn+89rpue+a/+vGS7TlXSURcAEPgILUEkLMSmn13bRR/9eohG9ExQtdPUnBW7NWLmSn26i466AIDARmgJQpe0jtTf7xqkF+8cqOSYCB04cUo/yf5COwr4uAgAELgILUGsrqPuNZfFqaLaqfvf2MBHRQCAgEVoCXIt7aGaObaf4qPs+rqwVH9atNXqkgAAaBJCSzMQ18qumWP6yTCkf3+2X0u25FtdEgAAbiO0NBP/c1mc7rn2UknSlLc26eDJUxZXBACAewgtzcivb+im1JTWKj5drUlzN6q6hmcWAQACB6GlGQkLsem52/qrlT1UX+w9oWc/3mV1SQAAuIzQ0sx0aNtC/3dLb0nSrI93at2eYxZXBACAawgtzdDN/S7RrQPby2lKD7zp0AkesggACACElmbq4e/1Upe4lsovOq0pb2+SafKARQCAfyO0NFMt7aF6dlx/hYfYtGxrof61bp/VJQEAcF6Elmas9yUxmjqyhyTp0UXbtC2/2OKKAABoHKGlmbv76k66rkc7VVY7df8bG2nzDwDwW4SWZs4wDP31h33VLsquXYdL9cf3afMPAPBPhBaobSu7Zo6tbfP/xuf7tWgTbf4BAP6H0AJJ0tVd4/SLIbVt/qe+s0kHTpRbXBEAAA2FWl0A/McDI7pp7Z5j2rj/pMa+sE4psZEuH9suKkK/G9lDl7R2/RgAANxhmEHSoKO4uFgxMTEqKipSdHS01eUErLzj5brp2f+q5HS128fGtbLrH1mDlJrS2vOFAQCCkjvv34QWnGXfsTJtPljk8ninKf3tk13aXlCiiDCbnh7bTzf2TvJihQCAYEFoIbT4XGlFte57fYNW7DgiSZo6sofuubaLDMOwuDIAgD9z5/2bibjwiFb2UL101yDdld5RkvT4B9s17Z3NqqpxWlwZACBYEFrgMaEhNv3x5t6antlTNkOa+0We7n7lCxWdqrK6NABAECC0wOPuvrqz/n7XILUID9HqXUf1gzlrlHecJdQAgItDaIFXXH95gubdk67E6AjtOlyq0bM/Vc6+E1aXBQAIYIQWeE3vS2K0cMLV6pUcrWNllRr393V678tDVpcFAAhQhBZ4VWJMhObdk67hl3/zUMbnlu9UkCxaAwD4EEue4RM1TlOPLd6mf6zOlSTd0v8SXdejncvHh4UYSr80TjGRYd4qEQBgAfq0EFr81mvr9unh/3ylGqf7f+3aRdn1lx/21dDurocdAIB/I7QQWvzaf3ce0T9W5+p0VY3Lx+QdP6WDJ09Jku5I66Dfj7pcLcJ5dBYABDpCC6El6JyuqtHjH2xX9pq9kqRObVvoyTH9NLBjG2sLAwBcFDriIuhEhIXo4e/10r/HpykpJkJ7j5Xr1ufX6K9Lt6uymq67ANAcuB1aVq1apczMTCUnJ8swDC1cuPC84/Pz83X77berW7dustlsmjRp0jnHzZ8/Xz169FBERIT69OmjxYsXu1samoGru8ZpyaRr9f3+l8hpSrM/2a3Rsz/VjoISq0sDAHiZ26GlrKxMqampmj17tkvjKyoqFB8frwcffFCpqannHLNmzRqNGzdOP/3pT7Vx40aNHj1ao0eP1pYtW9wtD81ATGSYnhrbT3PuGKA2LcK0Nb9Ymc+t1ourdjdpgi8AIDBc1JwWwzC0YMECjR492qXxQ4cOVb9+/fT000832D927FiVlZXp/fffr9935ZVXql+/fnr++eddem3mtDRPh0tOa9rbm7V8+2FJ0uDOsXry1lSlxLawuDIAgCsCbk7L2rVrNXz48Ab7MjIytHbt2kaPqaioUHFxcYMNzU+7qAi9lDVIj3+/j1qGh+jz3OO68elVevOL/TSwA4Ag4xehpaCgQAkJCQ32JSQkqKCgoNFjZsyYoZiYmPotJSXF22XCTxmGodsGd9AHE6/VFZ3aqKyyRr97e7N+Ndehqhom6QJAsPCL0NIU06ZNU1FRUf2Wl5dndUmwWIe2LTT35+maNrKHwkIMvfflIU2cu5HgAgBBwi+6cyUmJqqwsLDBvsLCQiUmJjZ6jN1ul91u93ZpCDAhNkP3DLlU3RKidM9rOVq8uUCGHHrmtn4KDQnYjA4AkJ/caUlPT9fy5csb7Fu2bJnS09MtqgiBbliPdnr+zgEKCzG0aHO+Jr7pUDV3XAAgoLkdWkpLS+VwOORwOCRJubm5cjgc2r9/v6Taj23uuuuuBsfUjS8tLdWRI0fkcDi0devW+u9PnDhRS5Ys0ZNPPqnt27fr4Ycf1vr163XfffddxKmhubuuR4Lm3DGwNrhsytcD874kuABAAHN7yfOKFSs0bNiws/ZnZWUpOztbP/7xj7V3716tWLHim9/EMM4a37FjR+3du7f+6/nz5+vBBx/U3r17ddlll+kvf/mLbrrpJpfrYskzGrNsa6F++e8cVdWY+l5qsp4ak8pHRQDgJ3j2EKEF3/HhVwX65b83qNpp6uZ+yXpqTD+F2M4O0wAA3wq4Pi2At93QK1Gzbh+gUJuhdx2H9Jv5X9I9FwACDKEFzcaNvRM16/b+CrEZWrDxoH5LcAGAgEJoQbNyY+8kPTeuNri8s/Ggpry1ieACAAHCL/q0AL50U58kmab0q7kb9faGAzIM6S8/6Csbc1wAwK8RWtAsjeqbJFOmJs516K2cAzIk/ZngAgB+jY+H0Gz9b99kzRzbTzZDmp9zQH94d4vVJQEAzoPQgmbte6m1wcUwpH9/tl/7j5VbXRIAoBGEFjR7N/e7RFd0ipUkrfj6sMXVAAAaQ2gBJA3tHi9JWrHjiMWVAAAaQ2gBJA3r3k6StGb3UZ2uqrG4GgDAuRBaAEk9EqOUGB2h01VOfZZ73OpyAADnQGgBVPtQzyHd6j4iYl4LAPgjQgtwxrAezGsBAH9GaAHOuLprnEJthnKPlmnfsTKrywEAfAehBTgjKiJMgzq1kcTdFgDwR4QW4FuGnllF9AnzWgDA7xBagG+p69eydvcxlj4DgJ8htADf0j0hSkkxEaqodmrtnmNWlwMA+BZCC/AthmHU321ZybwWAPArhBbgO+rmtdCvBQD8C6EF+I6ru8YpLMTQ3mPlyj3K0mcA8BeEFuA7WtlDNajjmac+c7cFAPwGoQU4h7ruuJ8wrwUA/AahBTiHunkt6/Yc06lKlj4DgD8gtADncFm7VrqkdaQqq51ax9JnAPALhBbgHAzD0JDudR8RMa8FAPwBoQVoxNBu3zz12TRNi6sBABBagEbULX3ef7xce1j6DACWI7QAjWhpD9XgznVLn1lFBABWI7QA5zG0G91xAcBfEFqA86jr1/LZnuMqr6y2uBoAaN4ILcB5XBp/ZulzjVNrd7P0GQCsRGgBzsMwjPq7LcxrAQBrEVqAC6ib1/LJjsMsfQYACxFagAu4qmtbhYfYdODEKe0+wtJnALAKoQW4gBbhoUrrwlOfAcBqhBbABUO6Ma8FAKxGaAFcMKxH7byWz3OPq6yCpc8AYAVCC+CCLnEtlRJbu/R5DUufAcAShBbABYZh0B0XACxGaAFc9O1+LSx9BgDfI7QALkrvEqfwUJsOnjylXYdLrS4HAJodQgvgosjwEF3Zpa0kVhEBgBUILYAbhp5Z+vwJ81oAwOfcDi2rVq1SZmamkpOTZRiGFi5ceMFjVqxYoQEDBshut6tr167Kzs5u8P2HH35YhmE02Hr06OFuaYDXDe1eG1q+2HtcpSx9BgCfcju0lJWVKTU1VbNnz3ZpfG5urkaNGqVhw4bJ4XBo0qRJGj9+vJYuXdpgXK9evZSfn1+/rV692t3SAK/rHNdSHdu2UFWNqU93HbW6HABoVkLdPWDkyJEaOXKky+Off/55de7cWU8++aQk6fLLL9fq1as1c+ZMZWRkfFNIaKgSExPdLQfwqdqlz/H659p9WrHjiDJ68XcWAHzF63Na1q5dq+HDhzfYl5GRobVr1zbYt3PnTiUnJ6tLly664447tH//fm+XBjTJ0O61/VpW8tRnAPApr4eWgoICJSQkNNiXkJCg4uJinTp1SpKUlpam7OxsLVmyRHPmzFFubq6uueYalZSUNPq6FRUVKi4ubrABvnBll7ayh9p0qOi0HHknrS4HAJoNv1g9NHLkSN16663q27evMjIytHjxYp08eVLz5s1r9JgZM2YoJiamfktJSfFhxWjOIsNDdFOfJEnSY4u3cbcFAHzE66ElMTFRhYWFDfYVFhYqOjpakZGR5zymdevW6tatm3bt2tXo606bNk1FRUX1W15enkfrBs7ntxndFRFm0xd7T+j9TflWlwMAzYLXQ0t6erqWL1/eYN+yZcuUnp7e6DGlpaXavXu3kpKSGh1jt9sVHR3dYAN8Jbl1pH45tKskacbibTpVWWNxRQAQ/NwOLaWlpXI4HHI4HJJqlzQ7HI76ibPTpk3TXXfdVT/+3nvv1Z49ezRlyhRt375df/vb3zRv3jw98MAD9WN+85vfaOXKldq7d6/WrFmjW265RSEhIRo3btxFnh7gPT+/tosuaR2pQ0WnNWflbqvLAYCg53ZoWb9+vfr376/+/ftLkiZPnqz+/fvroYcekiTl5+c3WPnTuXNnLVq0SMuWLVNqaqqefPJJvfTSSw2WOx84cEDjxo1T9+7dNWbMGLVt21br1q1TfHz8xZ4f4DURYSH6/ajLJUkvrNytAyfKLa4IAIKbYQbJLMLi4mLFxMSoqKiIj4rgM6Zpatzf12ndnuMa1SdJs+8YYHVJABBQ3Hn/9ovVQ0CgMgxD0zN7yWZIizbna+3uY1aXBABBi9ACXKTLk6J1e1oHSdIj732l6hqnxRUBQHAitAAe8OsR3RUTGabtBSV64wuW3wOANxBaAA9o0zJck0d0kyQ99eEOnSyvtLgiAAg+hBbAQ+5I66BuCa10orxKT3+00+pyACDoEFoADwkNsWl6Zi9J0mvr9mlHQePPzgIAuI/QAnjQ1V3jlNErQTVOU398/yueSwQAHkRoATzs9zf1VHioTZ/uOqYPtxZe+AAAgEsILYCHdWjbQj+7prMk6U+Ltup0Fc8lAgBPILQAXvDLoV2VEG1X3vFT+sfqXKvLAYCgQGgBvKClPVTTRtY+l2j2J7tUUHTa4ooAIPARWgAvublfsgZ2bKPyyhr9ecl2q8sBgIBHaAG8pPa5RD1lGNKCjQeVs++E1SUBQEALtboAIJj1bd9atw5sr3nrD+iR977S49/v6/KxLcJD1CmupRerA4DAYphB0kjCnUdbA750pKRC1z2xQiUV1W4fO3ZQih7/QR8ZhuGFygDAeu68f3OnBfCy+Ci7Hv5eLz354Q5VOV3/N8Kx0gq9uT5PHeNa6JdDu3qxQgAIDNxpAfzUa2v36g/vfiVJev5HA3Rj7ySLKwIAz3Pn/ZuJuICfujO9k7LSO0qSHnjzS205WGRxRQBgLUIL4Mf+8L89dW23eJ2qqtFP//mFCovp9wKg+SK0AH4sNMSmWbf3V9d2rVRYXKHx/1yvU5U8FgBA80RoAfxcdESYXs66Qm1ahGnzwSL9er5DTjcm9AJAsCC0AAGgQ9sWeuHOQQoLMbR4c4FmfvS11SUBgM8RWoAAMbhzrGacaU733Me7tHDjQYsrAgDfIrQAAeSHA9vr3iGXSpKmvLVJOfuOW1wRAPgOoQUIMFMyuuuGngmqrHHq56/mKO94udUlAYBPEFqAAGOzGXr6tn7qlRytY2WVGv/P9So5XWV1WQDgdYQWIAC1CA/VS1mD1C7Krh2FJZo416EaVhQBCHKEFiBAJcVE6u93DZI91KaPtx/WY4u3WV0SAHgVoQUIYKkprfXUmH6SpH+sztWra/daWg8AeBOhBQhwo/om6dcjukmSHnr3Kz27fKeC5DmoANAAoQUIAvdd11W/HFq7FPqpZV/r/y3YrOoap8VVAYBnEVqAIGAYhqbc2EOP3txLhiG98Xme7nktR+WV1VaXBgAeQ2gBgsid6Z00546BsofatHz7Yd3+9890rLTC6rIAwCMILUCQubF3ol7/WZpatwiTI++kfjBnjfYdK7O6LAC4aIQWIAgN7Birt39xldq3idTeY+X6/t/W6Mu8k1aXBQAXhdACBKlL41vpnV9eVd8597YX1+mT7YetLgsAmozQAgSxdlERevOedF1zWZxOVdVo/KvrNe+LPKvLAoAmIbQAQa6VPVT/yLpC3x9wiWqcpqa8vUnPfEQvFwCBh9ACNAPhoTY9eWuqJgyr7eUy8yN6uQAIPIQWoJkwDEO/zeihR0f3lu1ML5d7/5Wjiuoaq0sDAJcQWoBm5s4rO+r5H9X2cvlo22Hd//pGVXHHBUAAILQAzdANvRL1UtYghYfa9OHWQk2e96VqnMxxAeDfCC1AM3XNZfGac8cAhYUYeu/LQ5ry1iY5CS4A/BihBWjGrr88Qc+N668Qm6G3NxzQg+9uYVURAL9FaAGauRt7J+mpMakyDOn1z/brj+9vJbgA8Etuh5ZVq1YpMzNTycnJMgxDCxcuvOAxK1as0IABA2S329W1a1dlZ2efNWb27Nnq1KmTIiIilJaWps8//9zd0gA00c39LtGff9BXkvTKp3v1l6U7CC4A/I7boaWsrEypqamaPXu2S+Nzc3M1atQoDRs2TA6HQ5MmTdL48eO1dOnS+jFvvvmmJk+erOnTp2vDhg1KTU1VRkaGDh+m5TjgK2MGpejR0b0lSXNW7Nazy3dZXBEANGSYF/HPKcMwtGDBAo0ePbrRMb/73e+0aNEibdmypX7fbbfdppMnT2rJkiWSpLS0NF1xxRWaNWuWJMnpdColJUX333+/pk6d6lItxcXFiomJUVFRkaKjo5t6SkCz99J/9+hPi7ZJkqaO7KF7h1xqcUUAgpk7799en9Oydu1aDR8+vMG+jIwMrV27VpJUWVmpnJycBmNsNpuGDx9eP+ZcKioqVFxc3GADcPHGX9NFv83oLkl6/IPteuXTXIsrAoBaXg8tBQUFSkhIaLAvISFBxcXFOnXqlI4ePaqamppzjikoKGj0dWfMmKGYmJj6LSUlxSv1A83RhGFd9avrukqSHnlvq17/bL/FFQFAAK8emjZtmoqKiuq3vDyeXAt40gMjuunn13aRJP1+4Wa9s+GAxRUBaO5Cvf0bJCYmqrCwsMG+wsJCRUdHKzIyUiEhIQoJCTnnmMTExEZf1263y263e6VmALVz1qaN7KGKqhr9c+0+/Wb+lwqxGbq53yVWlwagmfL6nZb09HQtX768wb5ly5YpPT1dkhQeHq6BAwc2GON0OrV8+fL6MQCsYRiGpmf20m1XpMhpShPnOnTDzJWa9fFO7T9WbnV5AJoZt++0lJaWateub5ZC5ubmyuFwKDY2Vh06dNC0adN08OBBvfrqq5Kke++9V7NmzdKUKVP0k5/8RB9//LHmzZunRYsW1b/G5MmTlZWVpUGDBmnw4MF6+umnVVZWprvvvtsDpwjgYthshh67pY8iwkL0+mf79XVhqZ748Gs98eHXSk1pre+lJiuzb5LaRUdYXSqAIOf2kucVK1Zo2LBhZ+3PyspSdna2fvzjH2vv3r1asWJFg2MeeOABbd26Ve3bt9cf/vAH/fjHP25w/KxZs/TXv/5VBQUF6tevn5599lmlpaW5XBdLngHvKzpVpaVfFei9Lw/p011HVfeoIsOQ0ru01fdSkzWyd5JiWoRZWyiAgOHO+/dF9WnxJ4QWwLeOlFRo0aZD+s+Xh7Rh/8n6/WEhhoZ0i1dmarJG9ExQi3CvT51DgDFNU4ZhWF0G/AShhdAC+FTe8XK9t+mQ/uM4pO0FJU1+nVCboSu7tNVNfZKU0StBbVs1v8n2pmnqVFWNSiuqVVZRo9LT1SqtqD7zdbXKKqt1MQ/j/nZU+HZuML71HXfzRHWNU6UVNSr7Tp2lFTUqr9tXeeZ8KqpVWe2UYdRe71CbTaE2QyEhRu2vdftC6v7bUIjNJncjTnRkqOKjIhTfyq74KLviWoUrPspev7VtaVeIjeDkDwgthBbAMjsLS/SfL2vvwOy7iMm6ITZDV3aJPRNgEhUXBAGmusapfcfLtetwaf2271iZSk7XvtGXnHnDv5hQAtfYDCm2ZbjizoSalm7eEQwLtSm2RZhatwhXbMtwtWkZrjYtwtSmRe1/x7YIV2R4iJeqDy6EFkILYDnTNHWyvEpON/4Xc6K8Sh9uLdDizfnacvCbLtc2Q0rr3FY39U3Sjb0SFR/l3wHmdFWN9hwp087DJdp9uFS7jtQGlNyjZaqqce3PwzCkVuGhamkPVauI2l+j7KFqER7S5DsEdZfClHmOfQ2/dkeozait0x6ilvYzNdf/+s2+luGhamkPUWRYiJymVOM0VVXjVI3TVLXTPPOrU9U1Z3/tDqdpquhUlY6UVOhoaaWOlFToSGlF7a8lFTpeVuGTYBgRZqsNMS3ClRgToWsvi9Pwnglq36aF93/zAEJoIbQAAW/fsTIt3lwbYDYfLKrfbzOkwZ1jNapPkjJ6J6pdlHurlqprnCo+Xa2T5ZU6eapKRaeqVFRepZPllSo6Va2Tpyprvz7zvdNVNW69fsnpauWdKG/0zT8yLESXtmupy9pFqWu7Vuoc11KtI8Pqg0mrM1tkWIhsfHzhFTVOU8fLGoaZU25e54qqGp0or9SJ8iqdKKus/e+yKh0vr9TJ8srzhtOeSdEa3jNBN/RMUK/k6GY/v4fQQmgBgkre8XIt3pyvxZvz9eWBbwKMYUj2UNfbTZmmVFHt9EaJZ4mJDNNl7Vqp65nt0natdFm7VkqOiSSMBDnTNFVaUa0TZVU6UV6p4+WV+rqgRMu3Hdb6fccb3OVJionQ8MsTNKJngq7s0lbhbvx9DhaEFkILELTyjpfrgy35WrS5QF/mnWzy60TZQxXTIkytW4QpJjJMrSPDa7+OPPN1izDFRIYrIszm1r+EI0Jt6hLfSnGtwpv9v6BxtmOlFfp4+2F9tK1Qq74+2uAOTyt7qIZ0j9cNPRM0tHs7xUQ2j9YBhBZCC9AsHC+rVHlltVvHtAgPVXREqEJDmt+/aOFfTlfV6NNdR/XRtkIt23pYR0srmvxaoTZDmanJmp7ZU61bhHuwSu8jtBBaAAABxOk05ThwUh9tLdSyrYXaebi0Sa+TEG3XX3+Yqmu7xXu4Qu8htBBaAAAB7EKTeb8r92iZpr69SXuOlkmS7krvqGkjLw+IZdeEFkILAKCZOVVZoxkfbNOra/dJkrrEtdTMsf2UmtLa2sIuwJ33bz7UBQAgCESGh+iPN/fWqz8ZrIRou/YcLdP356zRzGVfq6rGN6vmvI3QAgBAELm2W7yWTrpWmanJqnGaemb5Tv1gzhrtauI8GX9CaAEAIMi0bhGu58b117Pj+is6IlSbDhRp1LP/VfanuXIG8HMiCC0AAASp76Um68MHhuiay+JUUe3Uw+9t1V0vf678olNWl9YkTMQFACDImaap19bt02OLt+l0lVPREaEa3f8St59jNXlEN0VFeLbpHauHCC0AAJxl95FSTX7T0eBxGO74/PfXu/28rwtx5/3bvWdxAwCAgHVpfCu9/Yur9FbOAeWdKHf7+Bbh1sYGQgsAAM1IaIhNtw3uYHUZTcJEXAAAEBAILQAAICAQWgAAQEAgtAAAgIBAaAEAAAGB0AIAAAICoQUAAAQEQgsAAAgIhBYAABAQCC0AACAgEFoAAEBAILQAAICAQGgBAAABIWie8myapiSpuLjY4koAAICr6t63697HzydoQktJSYkkKSUlxeJKAACAu0pKShQTE3PeMYbpSrQJAE6nU4cOHVJUVJQMw7C6HK8qLi5WSkqK8vLyFB0dbXU5PsE5c87BqLmdr8Q5c85nM01TJSUlSk5Ols12/lkrQXOnxWazqX379laX4VPR0dHN5gegDufcPDS3c25u5ytxzs2Fq+d8oTssdZiICwAAAgKhBQAABARCSwCy2+2aPn267Ha71aX4DOfcPDS3c25u5ytxzs2Ft845aCbiAgCA4MadFgAAEBAILQAAICAQWgAAQEAgtAAAgIBAaPFTs2fPVqdOnRQREaG0tDR9/vnnjY7Nzs6WYRgNtoiICB9We3FWrVqlzMxMJScnyzAMLVy48ILHrFixQgMGDJDdblfXrl2VnZ3t9To9yd1zXrFixVnX2DAMFRQU+KZgD5gxY4auuOIKRUVFqV27dho9erR27NhxwePmz5+vHj16KCIiQn369NHixYt9UK1nNOWcA/nnec6cOerbt299Q7H09HR98MEH5z0mkK+v5P45B/L1bczjjz8uwzA0adKk847zxLUmtPihN998U5MnT9b06dO1YcMGpaamKiMjQ4cPH270mOjoaOXn59dv+/bt82HFF6esrEypqamaPXu2S+Nzc3M1atQoDRs2TA6HQ5MmTdL48eO1dOlSL1fqOe6ec50dO3Y0uM7t2rXzUoWet3LlSk2YMEHr1q3TsmXLVFVVpRtuuEFlZWWNHrNmzRqNGzdOP/3pT7Vx40aNHj1ao0eP1pYtW3xYedM15ZylwP15bt++vR5//HHl5ORo/fr1uu6663TzzTfrq6++Ouf4QL++kvvnLAXu9T2XL774Qi+88IL69u173nEeu9Ym/M7gwYPNCRMm1H9dU1NjJicnmzNmzDjn+FdeecWMiYnxUXXeJclcsGDBecdMmTLF7NWrV4N9Y8eONTMyMrxYmfe4cs6ffPKJKck8ceKET2ryhcOHD5uSzJUrVzY6ZsyYMeaoUaMa7EtLSzPvueceb5fnFa6cczD9PJumabZp08Z86aWXzvm9YLu+dc53zsF0fUtKSszLLrvMXLZsmTlkyBBz4sSJjY711LXmToufqaysVE5OjoYPH16/z2azafjw4Vq7dm2jx5WWlqpjx45KSUm5YMoPdGvXrm3w5yNJGRkZ5/3zCRb9+vVTUlKSRowYoU8//dTqci5KUVGRJCk2NrbRMcF2rV05Zyk4fp5ramo0d+5clZWVKT09/Zxjgu36unLOUnBcX0maMGGCRo0addY1PBdPXWtCi585evSoampqlJCQ0GB/QkJCo/MXunfvrpdfflnvvvuu/vWvf8npdOqqq67SgQMHfFGyzxUUFJzzz6e4uFinTp2yqCrvSkpK0vPPP6+3335bb7/9tlJSUjR06FBt2LDB6tKaxOl0atKkSbr66qvVu3fvRsc1dq0DaS5PHVfPOdB/njdv3qxWrVrJbrfr3nvv1YIFC9SzZ89zjg2W6+vOOQf69a0zd+5cbdiwQTNmzHBpvKeuddA85bk5S09Pb5Dqr7rqKl1++eV64YUX9Oijj1pYGTyle/fu6t69e/3XV111lXbv3q2ZM2fqtddes7CyppkwYYK2bNmi1atXW12Kz7h6zoH+89y9e3c5HA4VFRXprbfeUlZWllauXNnom3gwcOecA/36SlJeXp4mTpyoZcuW+XwSMaHFz8TFxSkkJESFhYUN9hcWFioxMdGl1wgLC1P//v21a9cub5RoucTExHP++URHRysyMtKiqnxv8ODBAfmmf9999+n999/XqlWr1L59+/OObexau/qz4C/cOefvCrSf5/DwcHXt2lWSNHDgQH3xxRd65pln9MILL5w1Nliurzvn/F2Bdn0lKScnR4cPH9aAAQPq99XU1GjVqlWaNWuWKioqFBIS0uAYT11rPh7yM+Hh4Ro4cKCWL19ev8/pdGr58uXn/Yz022pqarR582YlJSV5q0xLpaenN/jzkaRly5a5/OcTLBwOR0BdY9M0dd9992nBggX6+OOP1blz5wseE+jXuinn/F2B/vPsdDpVUVFxzu8F+vVtzPnO+bsC8fpef/312rx5sxwOR/02aNAg3XHHHXI4HGcFFsmD19r9+cLwtrlz55p2u93Mzs42t27dav785z83W7dubRYUFJimaZp33nmnOXXq1PrxjzzyiLl06VJz9+7dZk5OjnnbbbeZERER5ldffWXVKbilpKTE3Lhxo7lx40ZTkvnUU0+ZGzduNPft22eapmlOnTrVvPPOO+vH79mzx2zRooX529/+1ty2bZs5e/ZsMyQkxFyyZIlVp+A2d8955syZ5sKFC82dO3eamzdvNidOnGjabDbzo48+suoU3PaLX/zCjImJMVesWGHm5+fXb+Xl5fVjvvt3+9NPPzVDQ0PNJ554wty2bZs5ffp0MywszNy8ebMVp+C2ppxzIP88T5061Vy5cqWZm5trbtq0yZw6dappGIb54YcfmqYZfNfXNN0/50C+vufz3dVD3rrWhBY/9dxzz5kdOnQww8PDzcGDB5vr1q2r/96QIUPMrKys+q8nTZpUPzYhIcG86aabzA0bNlhQddPULef97lZ3jllZWeaQIUPOOqZfv35meHi42aVLF/OVV17xed0Xw91z/vOf/2xeeumlZkREhBkbG2sOHTrU/Pjjj60pvonOdb6SGly77/7dNk3TnDdvntmtWzczPDzc7NWrl7lo0SLfFn4RmnLOgfzz/JOf/MTs2LGjGR4ebsbHx5vXX399/Zu3aQbf9TVN9885kK/v+Xw3tHjrWhumaZru3ZsBAADwPea0AACAgEBoAQAAAYHQAgAAAgKhBQAABARCCwAACAiEFgAAEBAILQAAICAQWgAAQEAgtAAAgIBAaAEAAAGB0AIAAAICoQUAAASE/w99xKSsEXT8WgAAAABJRU5ErkJggg==\n"
          },
          "metadata": {}
        }
      ],
      "source": [
        "import matplotlib.pyplot as plt\n",
        "plt.plot(thresholds, odd_ratios)\n",
        "plt.show()"
      ]
    }
  ],
  "metadata": {
    "colab": {
      "name": "mbp_aki.ipynb",
      "provenance": [],
      "include_colab_link": true
    },
    "kernelspec": {
      "display_name": "Python 3 (ipykernel)",
      "language": "python",
      "name": "python3"
    },
    "language_info": {
      "codemirror_mode": {
        "name": "ipython",
        "version": 3
      },
      "file_extension": ".py",
      "mimetype": "text/x-python",
      "name": "python",
      "nbconvert_exporter": "python",
      "pygments_lexer": "ipython3",
      "version": "3.9.12"
    }
  },
  "nbformat": 4,
  "nbformat_minor": 0
}